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ABSTRACT 


POLARIZATION OF THE LIGHT FROM 
THE 3 1 P-2 1 S TRANSITION IN PROTON BEAM- 
EXCITED HELIUM 

by 

MARTIN S. WEINHOUS 

Measurements of the polarization of the light from 
the 3 1 P-2 1 S (X 5016 A) transition in proton beam excited 
Helium have shown both a proton beam energy and Helium 
target gas pressure dependence. Results for the linear 
polarization fraction (at right angles to the proton beam 
and at .2 mtorr He target pressure) range from +2.6% at 
100 keV proton energy to -5.5% at 450 keV. The zero cross 
over occurs at approximately 225 keV. This is in good 
agreement with other experimental work in the field, but 
in poor agreement with theoretical predictions. The 
other experimental workers have used .2 mtorr as their 
lowest He target gas pressure, while in this work measure- 
ments have been made at He target gas pressures as low 
as .01 mtorr. The results have shown that the linear 
polarization fraction is still pressure dependent at .01 


mtorr. 



Me also have found a pressure dependence of photons 
per proton per He target _atom. We then conclude that 
experimental determinations of the linear polarization 
fraction have not yet been made under conditions which 
allow for strict comparison with theoretical predictions. 



SECTION I 


INTRODUCTION 

1*1 Purpose of the Investigation 

In the early 1900's it was known that spectral 
lines could exhibit polarization when magnetic fields 
were applied to the light source (the r and o components 
in a Zeeman effect spectrum) . However in about 1920 
the yellow mercury lines, 5770 A and 5791 A created in a 
gas discharge tube, were found to be weakly polarized 
even in the absence of a magnetic field. The light was 
polarized such that the maximum electric vector was par- 
allel to the path of current in the discharge tube. This 
fact led the investigators of the period., notably Skinner 
(1926) to test the hypothesis that the polarization was 
caused by the electron "beam” within the discharge tube. 
The hypothesis proved to be correct and so began the 
study of coll i s ionally- produced polarized atomic line 
radiation. Many such studies have been completed since 
that time, Pjg rci va 1 (1958) has published an extensive 
article on electron excitation of polarized atomic line 
radiation. Included in this article are Born approxima- 
tion calculations for excitation of the various magnetic 
substates of the target atom. More recently, 
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investigators have become interested in the polariza- 
tion due to proton impact. Bell (1961) has calculated 
theoretical excitation cross -sect ions for proton- 
excited helium. He has used both the Born and Distorted 
Wave Approximations for his calculations. The results 
of these calculations can easily be converted into an 
expected polarization of the light emitted after the 
collision. Two independent research groups have done 
just that. A Dutch group, Van Eck (1964) and Van Den Bos 
(1968) were able to compare these theoretical polarization 
values with experiments in the 5 to 35 keV and 1 to 150 
keV' proton energy ranges respectively. A second group 
working at the University of Giessen, Germany, has also 
done an experimental check of Bell’s work. Scharmann 
(1967), (1969) has also investigated the polarization of 
light emitted by He after proton impact in the energy 
range 100 to 1000 keV. Unfortunately, there is little 
overlap and only mediocre agreement between these two 
groups . 

It is the purpose of this thesis . then to make a 
detailed study of the polarization of the light emitted 
by He which has been excited by proton impact; to provide 
a confirmation of the work of either the Dutch or German 
research groups; and to extend the work in the direction 
of lower He target gas pressures in the hope of finding 
the "free atom" value for the polarization. Only then 
can one make a valid comparison with the theoretical 
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models which are available. 

1.2 Polarization Measurements 

The polarization of linearly polarized light is 
usually described by a quantity called the linear polari- 
zation fraction, denoted by the symbol it. Two quantities 
are required to calculate it, the light intensities with 
the electric field vectors parallel to and perpendicular 
to a "preferred" direction. The observation is made 
along a line which meets the "preferred" direction line 
at right angles. See Figure 1. Then 



obviously tt can range from -1 to +1. 

To measure' the polarization one must separate 
and measure the intensities of I f/ and Ij__. This can be 
done by a number of methods. In this work a glass 
laminated polaroid-type HN 32 sheet polarizer was placed 
on the observer's line of sight such that the light of 
interest passed normally through the polarizer and such 
that the polarizer could be rotated (about the line of 
sight) by 90°. If one now replaces the observer by an 
instrument capable of measuring intensities and then uses 
the polarizer to pass first E ;/ light and then Ej_ light 
for measurements of I ;/ and I respectively, one can then 
calculate ir. 




5 


In this experiment the "preferred” direction is 
the proton beam direction?, the light originates from a 
He gas filled target chamber, an interference filter 
selects the 5016 A line {of He) , and a photomultiplier 
measures the light intensity. See Figure 2. 

A detailed discussion of polarized light and its 
measurement is found in Appendix A. 

1.3 Astrophysical Interest 

Astronomy has its very roots in the observation 
of the heavens via visible light. As the science matured 
more and more information was extracted from that light. 
Just as the intensity, wavelength, and phase of the 
light. convey information to the observer, so does the 
polarization. In fact the magnetic field of some stars 
has been determined from measurements of the polarization 
of the Zeeman components of spectral lines. A similar, 
analysis of light from the sunspots on our own sun lead to 
a determination of the intensities and polarities of the 
magnetic fields associated with those spots. The radial 
polarization exhibited in the light from reflection 
nebula can be used to calculate the average particle 
size in the nebula. The correlation between the inter- 
stellar reddening of starlight and the amount of polari- 
zation of that light can be used to gain information about 
the magnetic field of our galaxy. The light from the 
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MEASUREMENT OF LINEAR POLARIZATION 



FIGURE 2 


t 



’’jet’ 1 of material emanating from the giant galaxy M87 
(in Virgo) is highly polarized and therefore is believed 
to be generated by synchrotron radiation. The astro- 
physicist thus needs to be informed as to what processes 
create polarized light and to what degree they create 
it. Only then can he successfully unfold his data. 

This work investigates one non-magnetic process which 
creates polarized light, a process which is certainly 
active in space, a proton collision with an atom. 

1.4 Previous Experimental Work 

Van Eck (1964) provides us with the first experi- 
mental determination of the polarization of the light 
resulting from proton impact on ground-state Helium. He 
used a Glann- Thompson prism as the polarization analyzer 
and a Leiss monochromator to isolate light from the 
transition being investigated. Since the intensity of 
reflected light from the grating in the monochromator de- 
pends upon the polarization of the light, Van Eck et al . 
had to take great care to separate the polarization of 
the light due to the atomic transition from that due to 
the instrumental polarization. His measurements covered 
a proton energy range of 5 to 35 keV (generated by a Von 
Ardenne ion source) and a Helium target gas pressure 
range of. .2 to 1 m torr . A sensitive McLeod gauge was 
used to determine the pressure of the Helium. 
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Van Den Bos (1968) essentially repeated Van Eck’s 
work {with some additions) . He extended the proton energy 
range up to 150 keV and added magnetic shielding to the 
collision chamber. The agreement between Van, Den Bos and 
Van Eck is very good (5 to 35 keV) . 

Scharmann (1967, 1969) used a very different de- 
tecting apparatus to find the polarization of the light 
from Helium which had been excited by proton impact. His 
detector used a sheet polarizer and interference filter 
rather than the Glan-Thompson prism and monochromator of 
the Dutch groups. He was also able to cover a very wide 
proton energy range of 100 to 1000 keV, The Helium tar- 
get gas pressure range in his study was .2 m torr to 5 m 
torr. 

1.5 Previous Theoretical Work 

Percival r s (1958) article on the "Polarization of 
Atomic Line Radiation Excited by Electron Impact" dis- 
cusses both the Oppenheimer-Penney Theory and Born approxi- 
mation methods for calculating the polarization. The 
Oppenheimer-Penney (O-P) theory is used to calculate the 
polarization of atomic line radiation when the cross- 
sect ons for exciting quantum states of the upper level 
are known. These upper level quantum states have a defi- 
nite orbital angular momentum component (the preferred 
direction being that of the electron beam). Percival 
pays particular attention to He and certain isotopes of 



Hg for which the nuclear spin is zero. In fact one can 
extract from his table of polarization formulae for the 
He multiplets (SL->-SL ' ) the expression 


IT (si 


- .a . & 


where the parameters G r h 0 and hj are determined by the 
values of S and L', then for our case (3 1 P^2 1 S), 

L' = 0; and from Percival's table G = 1, h Q = 1, h * 1.- 
Q 0 and Q x are of course the cross-sections for exciting 
the = 0, and = ±1 substates. This result is in- 
dependent of the exciting particle. 

Percival goes on to actually calculate via the 

Born approximation the Q for the excitation of the 

L 

3 J D states of He. His results however were in poor 
agreement with the experimental results available at the 
time. Percival questions the suitability of the approxi 
mate wave functions used rather than the validity of the 
Born approximation. 

Bell * s (1961) work uses both the Born and Distor 
ted Wave Approximation methods to calculate the cross- 
sections for the process 

H + He (l s l ) S — * H* + Hi (is Ip ) 1 P 

He used product wave functions for the ground. state of 
the target system and excited state wave functions 
portional to %( f >) * %( r ,) f, (O] 


pro- 



The actual forms of the contain "adjustable” param- 
eters which are chosen so as to obtain particular oscil- 
lator strengths. The calculations of the.-Qg * s are 
carried out by numerical methods on high speed computers 
Theoretical details are given in Appendix B. 



SECTION II 

EXPERIMENTAL APPARATUS 

2.1 Introductory Description of the 
Experimental Apparatus - 

In order to implement the polarization measure- 
ment discussed in section 1.2, apparatus was assembled 
which would provide a relatively stable proton beam, 
a confined region of Helium target atoms, and a method 
for measuring the intensity of the two linear polariza- 

o 

t ion components of the X 5016 A line of Helium. This 
apparatus was assembled in the physics department of 
the University of New Hampshire. The three main com- 
ponents of the experimental apparatus were a Van De Graaff 
positive ion accelerator and associated beam tube optics, 
differentially pumped target chamber and vacuum system; 
and the polarization detection system. The proton beam 
is of course produced by the accelerator, the Helium 
target atoms are isolated within the differentially 
pumped target chamber, and the polarization detection 
system will isolate and measure the intensity of the 
linear polarization components of the He line.. Each of 
these systems is described in detail in the following 
. paragraphs. 
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2-2 The Van De Graaff Accelerator 

Our accelerator is a model PN-400 manufactured 
by High Voltage Engineering of Burlington, Massachusetts. 
It is capable of producing positive ion beams within the 
energy range of ^100 to ^450 keV. In our work we confine 
ourselves to proton beams of ^1 to ^19 u A 1 s current. 
Within the accelerator, the protons are generated in a 
radio-frequency source bottle. Hydrogen gas is continu- 
ously leaked into the source bottle through a palladium 
leak (maintaining very high purity) , where it is ionized 
by a radio-frequency discharge. This source bottle is 
located at the high potential end of the accelerator tube. 
A small canal (beryllium) connects the source bottle and 
accelerator tube. The protons then enter the accelerator 
tube at a rate dependent upon both the hydrogen gas 
pressure in the source bottle, and the magnitude of the 
positive probe voltage within the bottle. Once into the 
accelerator tube, the protons are confronted with a 
focusing electric field and then an accelerating electric 
field. These fields are maintained by a focus plane 
and a series of equipotential planes. The voltages with- 
in the Van de Graaff are due to the relocation of charges 
by the motor driven belt. So the potential difference 
seen by the protons depends upon the amount of charge 
carried by the belt, and is adjustable. The proton beam 
energy then depends only on the potential difference 
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through which the proton accelerates. 

In -fact,- not only protons appear in the- beam. 

Any positive ions created in the source bottle will be 
accelerated in the beam. One will therefore get at 
least K , / and * 'Any impurities found in the 

bottle may ionize and produce accelerated positive ions. 
Rough measurements have shown our H + yield to be approxi- 
mately 10-40 per cent of total beam current. Since this 
beam is later magnetically analyzed, its content is not 
important so long as sufficient h"*" is present. Diagrams 
of the accelerator are shown in figures 3 and 4. 

2.3 The Accelerator System 

A large number of accessory systems are required 
by the Van de Graaff. The accelerator tube and beam tubes 
must be maintained at low pressures, and the ion beam 
must be steered and focused as well as energy stabilized. 
Our experiment also requires accessory systems to dif- 
ferentially pump the target chamber, and to accurately 
measure, the pressure within the target chamber. The ac- 
celerator is also used for another research project as 
well as for teaching. Some of our equipment has there- 
fore been designed around these other requirements. 

The vacuum system and beam tube arrangement is 
shown in fig. 5. Both the main pump and left port pump 
No. 2 are NRC six inch diffusion pumps backed by Welch 
mechanical pumps. Typical operating pressures are; at 
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BEAM OPTICS OF THE VAN DEGRAAFF ACCELERATOR 
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the main pump ^6-8x10" 6 torr; and at the No. 2 pump 
^10x10“ 6 torr. Left port pumps No. 1# and No. 3 are two 
inch diffusion pumps again backed by mechanical pumps. 

The center port beam tube has no pumps of its own, how- 
ever, it is short enough not to need one. The only time 
we use the center port is during the initial tune up of 
the Van de Graaff. The right port beam tube (details not 
shown in fig. 5) is used generally for teaching experi- 
ments and not for research, and will hence be ignored. 

A magnetic analyzer is used to select for a mono- 
energetic proton beam.. It is well known that when a 
charged particle moves with its velocity perpendicular 
to a constant homogenous magnetic field, the particle 
will follow a circular path while it remains within the 
magnetic field. The Lorentz force on such a particle 
is 

P . 8 2.3. 

where in our case z = 1 and q = e. The force equation 
is then 

I - v* s 

where for our nonrelativistic case 

J 



using 2.3.3 in 2.3.2 one arrives at 


„ my" 
r 


2.3.2 


2.3.3 
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- --fv -% {'iEnY'-. - - 2.3.4 

<?8 ^ 

This last equation shows that for a constant magnetic 
field B, and a constant beam energy E , the radius of 
curvature for any ion depends upon the square root of 
the ion’s mass. The magnetic analyzer then can bend the 
different ions in the beam into different paths and thus 
isolate them. In practice, in our system, the accelerator 
will be tuned to a specific energy and beam current in 
the center port. The analyzing magnet will then be turned 
on with a current known to be insufficient to bend pro- 
tons (the most easily bent positive ion) into the left 
port. The magnetic current will then be gently increased 
until our various indicators show a beam in the left 
port tube. The current supplied to the analyzing magnet 
is regulated by an Atomic Laboratories Inc. power supply 
and regulator (Model C) . The regulator holds the current 
(and hence magnetic field) constant to one part in 10 5 .- 
An Energy Stabilization System is used to prevent 
changes in the ion beam energy. The energy of the ion 
beam, as it exits from the Van de Graaff accelerator, 
depends upon the difference in potential between the ion 
source bottle and ground. This potential difference in 
turn is controlled by the quantity of charge located at 
the high voltage terminal of the accelerator. High 
Voltage Engineering’s design of the PN-400 accelerator 
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includes a corona probe extending inward from the pre- 
sure tank wall toward the high voltage terminal. This 
corona probe is tipped with an array of needle-like points, 
which increase the probe ' s . efficiency in draining charge 
from the high voltage terminal to the pressure tank and, 
therefore, ground. It is then clear that the beam energy 
is a function of the rate of charge leakage through the 
corona probe. The High Voltage Engineering Inc. Corona 
Stabilizer takes advantage of the above. The corona 
probe is not connected directly to tank ground, but 
rather to the plate of a type 4-125A vacuum tube. One 
can then control the charge drain from the high voltage 
terminal by controlling the conduction of the 4-125A, 
the corona stabilizer does exactly that and hence controls 
the beam energy. 

When a beam has been magnetically analyzed, and 
directed down the left port, it encounters two probes 
(an insulated vertical slit) within the beam tube. If 
the beam energy drifts upward (due to drifts within the 
accelerator) , then the beam will be bent less by the 
magnet and it will impact upon the center (or High energy) 
probe more heavily than upon the outside (or low energy) 
probe. This imbalance in probe currents is detected by 
the High Voltage Engineering Corona Stabilizer Ampli- 
fier, and this amplifier in turn sends a negative feed- 
back signal to the (grid and cathode) 4-125A. This 
signal changes the tube's conductance, so. as to correct 
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the beam energy, i.e., to have equal impact on the 
probes (this centers the correct energy beam in the tube) . 

Energy calibrations of the Van cle Graaff accelera- 
tor are done with a High Voltage Engineering generating 
voltameter. This instrument is used to measure and dis- 
play (digitally) the potential of the high voltage 
terminal (and, therefore, the beam energy). The unit 
consists of a chopper (or rotor) and stator plate (both 
eight sectioned) . The unit is located within the Van 
de Graaff pressure tank near the high voltage terminal. 

As the chopper rotates, it alternately exposes and shields 
the stator from the high voltage terminal. The voltage 
induced on the stators is then a chopped D.C. or roughly 
triangular A.C. which is proportional to the high vol- 
tage terminal potential and, therefore, proportional to 
the beam energy. The rectified output of the generating 
voltameter is now connected to a digital voltameter for 
a fast, easy, and accurate readout. We estimate an 
instrumental accuracy of from 1 to 2 per cent. This, 
however, assumes a linear response, a "good" calibration, 
and a focus voltage setting which remains at its cali- 
bration value. 

In practice, the linearity has been verified for 
a two point calibration, and the focus voltages used 
do not vary more than ^5 kV. The calibration procedure 
for the generating voltameter involves the use of the 
F I9 (p,ay)0 1G nuclear resonance. The cross-section for 
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this reaction shows two distinct resonances, one at a 
(laboratory) proton energy of 340 keV P and the second at 
a proton energy of 484 keV. To perform the calibration 
experiment, a small aluminum wafer is exposed to con- 
centrated hydroflouric acid for 15 min. This provides 
us with a thin target for use at the end of the beam tube. 
A 2 inch Nal scintillation detector is used to detect, 
the 6 MeV y rays emitted from the fluorine. The pulses 
from the detector were Amplified and sent to a single 
channel analyzer and then to a scalar. The scalar is 
generally set to repeat 10 second data acquisition periods 
and ^5 second display periods. One then impacts a pro- 
ton beam on the fluorine target, adjusting the energy of 
the beam to a value less than that required for a reso- . 
nance. The "background" on the scalar is then noted. 

One then gradually increases the beam energy (by increas- 
ing the magnetic current, and belt charge) until a peak 
is reached, again the number of 7 - ray counts is noted. 

The resonance occurs at that point in energy where the 
number of counts on the scalar is just the average of the 
background and peak values. The Van de Graaff’s energy 
is adjusted so that t;.. scalar is showing just that 
average value, and the digital voltameter is set to the 
resonance valuer We will generally use the 340.5 KeV 
value for our calibration and the higher energy resonance 
for a linearity check. Using this calibration procedure, 
an overall energy readout accuracy of ^8 KeV appears 
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appropriate. 

An. Electrostatic Focusing System was installed 
as ah improvement for the pulsing system used by Dotchin 
et. al. for mean life studies.. It has, however, become 
useful in this work due to its ability to focus the. beam 
through the differential pumping slits (to be discussed 
later) . The system was designed and constructed by 
D. L. Keator as a course project. In operation, one 
simply monitors the beam current (by means of a Farady 
cup) and adjusts the voltages supplied to the electro- 
static focusing system so as to maximize that current. 

The differentially pumped target chamber, shown 
in figure 6, is used to provide the relatively high 
pressures of He (<2xl0“ 3 torr) required for the experi- 
ment without filling the beam tube with gas and adversely 
influencing the operation of the accelerator. To main- 
tain a pressure differential between the target chamber 
and beam tube, one continually leaks the target gas into 
the chamber which is then pumped through narrow slits 
into the beam tubes. The slit size is a compromise be- 
tween maximizing beam current, and minimizing the flow 
of gas into the beam tube. In our case, the upstream 
slits (nearer the accelerator) are also a part of the 
pulsing system used by Dotchin et. al. There are three 
slits in that group, each a horizontal opening %1 x ^10 mm 
separated by a distance of 'vlS mm. The downstream slit 
is adjustable and set such that it will not intercept any 
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of the beam (typically 10 mm 45 mm) , During operation, 
the highest pressure used in the target chamber is 
^1 . 5x10 3 torr , a simultaneous pressure . measurement 
'WO cm upstream from the triple slit will typically indi- 
cate a pressure of 6xlQ” 6 torr, a ratio of 250:1. The 
lowest pressure used in the target chamber is W0" 5 torr 
at which time the pressure in the beam tube is Wxl0“ 6 
torr giving a ratio of W:l. Again, differential pump- 
ing is required because the beam tube and accelerator 
must be kept at pressures <10“ 5 torr? and because we want 
to spatially localize the beam gas collisions. 

The target chamber is a brass cylinder with 
5 cm I.D. and 17.7 cm length. Two glass rectangular 
windows have been attached on opposite sides of the 
cylinder walls with epoxy cement. The chamber is 
oriented such that one may look horizontally through both 
windows. See again figure 6. The window sizes are 
both 13 cm x 2.3 cm. This target chamber is used both 
for the mean life studies of Dotchin et. al. and for 
our polarization studies, each group using one window 
of the chamber. There are two vacuum couplings on the 
top of the chamber, one connected through a series of 
regulators and valves to the gas supply, and the other 
through valves to our pressure gauge. See figure 7. 

The gas supply system is quite conventional with 
a two stage regulator reducing the tank pressure to a 
value of WO psig. The gas then flows through >a 
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flexible hose to another regulator. This second regu- 
lator is a Ma the son Co. vacuum type, capable of regula- 
ting its output from v5Q - 750 torr. We. typically run 
at ^300 torr. Next, the gas encounters a Hoke Micromite 
fine metering valve, which is used to control its flow 
into the evacuated target chamber. I would add that with- 
out the vacuum regulator, we would only be able to use 
the first turn of the 18 turn metering valve, whereas, 
with the vacuum regulator's pressure reduction, we are 
able to use ^1-1/2 turns of the valve to achieve the 
desired pressures within the target chamber. 

The pressure measuring system consists of a num- 
ber of vacuum pressure gauges spread about the beam tube 
and target chamber, see figures 5 and 7. The gauges 
located at the main pump, left port pump No. 1 and left 
port pump No. 2 are the cold cathode type usable from 
^1x10 6 torr to 150x10” 6 torr . The vacuum gauge on left 
port pump No. 3 is the ionization type usable from 
vlQ 9 to ^10 5 torr. The above gauges are used for general 
maintenance of the system and not for data acquisition. 

An extremely accurate and sensitive MKS Instruments Inc. 
Baratron capacitance manometer is connected to the 
target chamber for pressure data measurements. The 
manufacturer claims a general accuracy of ,1 per cent, 
and repeatability of .005 per cent. The method employed 
by the designers is to compare the unknown pressure 
^x a smaller known reference pressure p . <■ This 



comparison is done in our model 77 H-l pressure head. 

The head is divided by a highly stressed thin flat 
metallic diaphragm into two chambers, one for reference 
pressure, one for the unknown pressure. When the 
reference and unknown sides differ in pressure, the 
diaphragm is deformed, but the diaphragm is part of an 
A.C. capacitance bridge. Therefore, any imbalance .in 
pressure will be transformed into an imbalance in the A.C. 
bridge. Our type 77M-XR indicator translates this bridge 
imbalance into a pressure reading. One great advantage of 
this arrangement is the ability of the indicator to null 
(or balance) out (digitally) the gross pressure difference 
and to then use the full sensitivity of the instrument on 
the finer levels of imbalance. One can then achieve five 
significant digits in the readout. 

In order to achieve high accuracy, one must know 
the reference pressure , or reduce to a value low 
enough to cause negligible error. Recall the pressure 
indicated by the instrument is the difference in pressure 

P x ~ P r‘ We are usin 9 a Varian Vaclon pump (2l/sec.) to 
maintain a very low P^. We estimate from the current 
drawn by the Vaclon pump, P r to be <10“ 6 torr and would 
therefore have <10 per cent error in our lowest pressure 
data. 

2.4 The Polarization Detection System 

This system is used to detect the intensity of 
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both the vertically and horizontally linearly polarized 
components of the 5016 A (3^P - 2^S transition) light 
from 'the helium target gas. The system is shown in 
block diagram in figure 8 and in 1:1 scale in figure 9. 

The polarization analyzer consists of a type 
HN-32 polarizer which has been carefully mounted in a 
holder and aligned such that it can be rotated from one 
terminal position by 90 degrees to another terminal 
position and then back, etc. This arrangement allows 
the experimenter to set the analyzer (at one terminal 
position) to pass light whose electric field vector is 
vertical, and to then rotate the analyzer (blindly) to 
its other terminal position where it will pass only 
horizontally polarized light. These rotations are in' 
fact done by the experimenter in the darkened accelerator 
room. Thus the need for the two (90 degrees apart) 
terminal positions. 

Reflections from the surfaces of the polarizer 
are according to the manufacturer isotropic and <4 per 
cent; they can therefore be ignored. As previously 
stated, the ratio of transmissions for the desired: 
undesired components is ^1.5xl0 4 . This arrangement is 
then quite suitable for alternate analysis of the 
orthogonal components of linearly polarized light. 

o 

The 5016 A line of He is selected by an in- 
terference filter. This type of filter is a device 
which will transmit (by constructive interference) 
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only certain wavelengths of light. The transmission 
is polarization form independent, and will therefore 
not effect our measurement. This polarization form in- 
dependence was tested and confirmed by a rotation of 
the entire detector and a comparison of the polariza- 
tion results for the two detector positions. This test- 
ing and confirmation process is described in section IV. 
The filter used for this work was manufactured by Spectra 
Films Inc., Winchester, Mass. It has a peak transmittance 
of 60 per cent at 5018 A, and a full width at half maxi- 

o 

mum of 18.8 A for normal incidence. The slight difference 
in wavelength between our line {5016 A) and peak is not 
significant, it only reduces the percentage transmission 
for our, line to 58 per cent, a 2% loss. See figure 10. 

The need for such a filter becomes apparent when 
one examines the He spectrum. There are a large number' 
of prominent lines in the spectrum. Table 1 lists some 
of those lines near our 5016 A line. Also, a partial 
energy diagram for He is shown in figure 11, and figure 
12 shows a spectral scan in the region around 5016 A. The 
lines nearest 5018 A are at 5047 A and 4922 A. It is 
necessary that our interference filter not transmit these 
or any other lines. We have no problem when the Helium 
light is normally incident upon the filter, the 5047 A 
line is attenuated by ^99 per cent, and the 4922 A line 
is attenuated by >99 per cent. However, all the light- 
incident upon the filter is not normally incident, and 
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TABLE 1 


- 5016 A 


Wavelength (A) 

. Transition 

(5016 A 

4713 

4 3 S-2 3 P 

303 

4859 (He II) 

8-4 

157 

4922 

4 1 D-2 1 P 

94 

5047 

4 1 S-2 1 P 

-31 
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PARTIAL TERM DIAGRAM OF He 
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The a'P^S transition is indicated by a heavy line 

FIGURE 11 
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Xp 0ak is in fact/ a function of the angle cf incidence. 
According to. Baird-Atomic, a manufacturer of interference 
filters, the wavelength of peak transmittance is lowered 
as the angle of incidence increases, such that 



where 0 is the angle of incidence, A Q is the wavelength 
of peak transmittance for the angle of incidence, A n is 
the wavelength of peak transmittance for normal incidence 
and r) is the effective index of refraction of the filter. 
A short computer program was written to compute A 0 as 
a function of both 6 and rj* The program is listed in 
Appendix C . 1 , and the results are shown in Table 2. 

In our detector (see again figure 9) the largest 
possible incidence angle is V7°. Referring to Table 2, 
we then see that the worst possible case (n=l) has a 

o 

*peak ” 4980 A « . This worst case is still quite good 
since ^99 per cent attenuation is achieved at only 

o 

^22 A from A peak (see again Figure, 10) . Our interference 
filter is then quite sufficient to isolate the 5016 A 
line of helium. V 

A photomultiplier tube and associated electronics 
are used to measure light intensity. We refer again to 
Figures 8 and 9 where the tube, its housing and the 
various electronics are shown. The tube is an EMI 
6256S, ^2 inch diameter, S (Q) , S b C g photocathode, and 
very low (^2xl0~ 9 A) dark current at the operating 



TABLE 2 


PROPERTIES OF THE INTERFERENCE FILTER 


Index of Refraction 

1.0 

ii 

it 

1.4 

H 

«* 

1.6 
ii ■ 

ir 


Angle of Incidence (deg) 

o 

0 

4 

8 

0 

4 

8 

0 

4 

8 


X 9 (A) 

5018 

5005 

4969 

5018 

5011 

4993 

5018 

5013 

4998 
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voltage of -1850 volts.. The tube is wrapped with black 

tape and housed in an aluminum light-tight container 

with a single aperture in front of the photocathode. 

The tube base is wired as shown in Figure 13, and power 

is supplied by a Power Designs Pacific Inc., Model 2k-10 

high voltage power supply. The Spectral response of the 

0 

photocathode is near its peak at 5000 A (10 per cent 
quantum efficiency) and is well suited for measurements 

O' 

on the 5016 A line. The tube was usually operated at 
room temperature excepting very warm days when ice water 
cooling was necessary to bring background count rates 
back down to ’'normal.*’ 

Signal pulses from the tube were fed into a C.l. 
1416 Amplifier, these amplified pulses were next passed 
through a C.l. 1430 single channel analyzer which was 
used as a discriminator to reject pulses of less than a 
preset magnitude. The Scalar output of the single channel 
analyzer was passed to a C.l. 1470 scalar for counting. 

It was found that the. discriminator setting had little 
effect upon the final signal to noise ratio, therefore, 
a setting was chosen to give a reasonable counting rate. 
This same effect was noticed by Pegg (1970) , who used 
this same tube and similar electronics. 

A second photomultiplier tube was used briefly 
during this research, it was mounted directly on the 
second window of the target chamber and monitored the 
helium light output for normalization purposes. The 
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tube was an RCA 8575 run at -2500 volts, at room tem- 
perature'.. The arrangement is shown in figure 14. This 
system was used only long enough to confirm that light 
normalization and current normalization gave the same 
value for the linear polarization fraction, it was then 
abandoned so as to minimize the changeover time between 
the experiment of Dotchin et. al. and this work. 

Due to the fact that the proton beam from the 
accelerator is not stable, we must use a Faraday cup 
to determine just how much beam we have had in any counting 
period. The cup is really just an extension of the beam 
tube, of length 14.5 cm, with an oval entrance aperture 
(^10 mm x ^20 mm) . The downstream end is closed by a 
double end window of Vycor glass. See figure 15. We 
depend upon the length of the cup, and the entrance 
aperture to contain secondary electrons. A comparison of 
beam current measurements was made for two lengths of 
the cup with and without the aperture, since the results 
were essentially identical, we concluded that secondary 
electrons were not escaping and that the cup was contain- 
ing the proton beam. We therefore felt there was no 
need for a suppression grid. The electronics for the 
cup were shown in figure 8 , and consist of Keithly Model 
610 Microammeter , a homemade voltage to frequency con- 
vertor, a C.I. invertor, and a Mechtronics 700 scalar. 

In use, the microammeter is connected directly to the 
Faraday cup, the D.C. chart recorder output of the 
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LIGHT NORMALIZATION SYSTEM 
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microammeter is sent to the voltage to frequency conver- 
tor (built by L. W. -Dot-chin from plans in -the General 
•Electric Transistor Manual, General Electric Co., Syracuse, 
N. Y. 1964) . The pulses from the convertor are inverted 
to match the input requirements of the scalar where they 
are counted. The number of counts appearing on this 
scalar is directly proportional to the number of protons 
passing through the target chamber, and can therefore be 
used as a normalization base. 
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SECTION III 
. EXPERIMENTAL METHOD 
3,1 Introductory Description 

The methods used in this research have been 
evolved so as to maximize the usefulness and effective- 
ness of the apparatus with which we have worked. Methods 
were evolved to deal with the irregularities in the proton 
beam from the Van de Graaff accelerator. In fact two 
methods were developed and , used simultaneously iand gave 
the same value for the linear polarization fraction. The 
need for these normalization methods and their implementa- 
tion will be discussed in the next paragraph. 

The actual data taking process is dependent upon 
the type and settings of electronics used. One must re- 
main within the useable limits of the electronics while 
at the same time maximizing the signal to noise ratio. 

Our methods for achieving this end are described in 
paragraph 3.3. 

Finally the actual step by step data taking pro- 
cess is outlined in paragraph 3.4. 
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3.2 The Need for Normalization 

The Van de Graaff accelerator , being mechanically 
powered, is more susceptible to irregularities than most 
other pieces of apparatus. The proton beam from the 
accelerator typically undergoes three different types of 
transition, irregular oscillations of beam current, long 
term upward or downward drift of beam current, and 
sudden M./2 to 1 1/2 second complete cessations of beam. 

The irregular oscillations in the beam current 
are probably due to a number of factors including; 
irregular transport of charge by the belt within the 
accelerator? fluctuations in the probe voltage (the probe 
voltage supply is powered by a generator driven by the 
belt) ; and overcorrectioris in the negative feedback x 
loops of the energy stabilization system. In so far as 
our experiment is concerned, high frequency oscillations 
are of no importance, they would simply average over the 
>1 second observation period used in taking data. Low 
frequency oscillations must be accounted for in order to 
correctly compare intensity measurements taken at differ- 
ent times. 

Long term upward or downward drifts of beam 
current must be dealt with in the same manner as the low 
frequency irregular oscillations, especially if these 
drifts are noticeable within vL second periods. These 
long term drifts are most, likely due to changes in the 
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rate of flow of H 2 through the paladium leak. 

The sudden cessations of beam are usually due to 
sparks within the accelerator tube, sparks among the 
equipotential planes, or a spark from the high voltage 
terminal to the tank wall. Occasionally, the beam is lost 
for a short time when the energy stabilization system is 
unable to correct for an excess or lack of charge on the 
high voltage terminal. The beam energy will then be too 
high/low for the magnet to steer the beam into the left 
port beam tube and the beam will be lost until the charge 
situation corrects itself. In any case, whenever the 
beam is interrupted, the data taken at that time is not 
used, the machine is reset if necessary and new data is 
taken. The problem is to determine when an interruption 
has occurred. The problem is solved by making the experi- 
menter a part of the apparatus. While taking data, the 
experimenter stands at the end of the left port beam tube 
and stares at the Vycor end window of the Faraday cup. 

In striking the Vycor, the proton beam generates a dis- ' 
tinct blue light. With practice, the experimenter can 
spot beam interruptions (as an interruption of the blue 
light) which last for as short a time as >1/10 second. 

The data is then not used, the scalars are reset and 
new data is taken. 

Again, normalization is required so that one may 
compare intensity measurements made at different times 
even though the beam current is not constant in time. 



47 


Also, the experimenter is part of the apparatus not only 
in the sense of making adjustments, but also as a sensor- 
using' eyes and even ears to detect irregularities in the 
Van de Graaff performance. 

We have used two different normalization systems 
to monitor beam current. In the first case, we monitored 
the irregularities in the light output of the helium 
target gas. This system was based upon the assumption 
that the light output was directly proportional to the 
number of excitations which had occurred within the tar- 
get chamber. The second system was a bit more direct, 
we used a Faraday cup to collect the beam after its 
passage through the target chamber (assuming very little 
loss) and integrated the beam current. The, first system 
then normalizes to the number of excited helium atoms 
(assuming no saturation occurs) , while the second system 
normalized to the number of protons passing through the 
target chamber. We formulated the hypothesis that the 
number of excitations should be directly proportional to 
the number of protons passing through the chamber, i.e» 
that calculations based upon data taken simultaneously 
with both systems should give the same result. This 
was in fact the case demonstrating that our two nor- 
malizing systems were equivalent. Once we were satis- 
fied as to that equivalence, the light normalization 
system was deleted from the experiment. Its continued 
use would have created large scale problems durino the 



changeovers between this work and that of Dotchin et al .. 
All the data and results presented in this work are 
based upon the current normalization system. 

3„3 Setting the Detector System 

The system is "set" by properly adjusting all 
discriminators, amplifier gains, and miscellaneous other 
parameters.* 

A very simple and direct method is used to time 
the runs, in fact, all runs have the same duration. As 
shown in Figure 8 , an Ortec 48 0 pulser drives a pre~ 
se table Mectronics 7 02 scalar/timer. Pulses are gen- 
erated by the pulser in syncronization with the A.C. line 
frequency, i.e. 60 pulses/second. The Mectronics scalar/ 
timer is preset to turn off at 80 counts, (or 80/60 of a 
second = 1.3 second) . When the scalar/timer shuts off, 
it will sinzultaneously (i.e. in a time which is short 
compared with the average time between data pulses to 
the other scalers) shut off, via built-in gating cir- 
cuits, the other two scalars. The ' counting time of 1.3* 
second was in part forced upon us by the data count rates 
and by the tendency of the Van de Graaff to spark quite 
often. We wanted a counting time of the order of 1 to 
1 1/2 seconds in order to get convenient count rates and 
yet still be short enough so that the probability of a 
spark (beam interruption) would be small, 1.3 second 
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was the best available compromise . 

The amplifier gains, discriminators, and photo- 
multiplier tube high voltage have been set for single 
photon counting and zero dead time correction. Typical 
count rates are shown in table 3, 

The scalars used are capable of count rates on 
the order of 10 4 -10 5 counts/second before dead time cor- 
rections are needed. 

The signal count rate for the normalization 
scalar is controlled by the range switch on the Keithly 
microammeter . The output of the Keithly is the same per- 
centage of 10 volts as is the indicator of full scale. 

We limit, by choice of range, the percentage of full 
scale such that the input voltage to the voltage to fre- 
quency converted does not exceed 2 volts, which in turn 
limits the count rate. 

The criteria for the polarized light intensity 
scalar is that the tube be operating in the single photon 
counting mode and that extraneous noise not be counted. 
The count rate is then controlled by tube voltage, am- 
plifier gain, and the single channel analyzer's dis- 
criminator. 

The sources of background counts are noise pulses 
(both thermal and light noise) in the case of the polar- 
ized light intensity scalar; and zero offset of the 
Keithly microammeter in the case of the normalization 
scalar. To minimize background counts in the first case, 



Scalar 

Scalar/timer 
Light intensity 
Normalization 


TABLE 3 

SCALAR COUNTING RATES 
Signal + 

Background rate 

60/sec 

100 to 5000/sec 
M.3 00/sec 


Background rate 
0 

10 tg 50/sec 
^7/sec 



one darkens the accelerator room, uses lead bricks to 
shield the tube from tank x-rays, and when necessary, 
cools' the tube. In the case of . the normalization scalar 
one can "adjust" the background count rate by adjusting 
the zero of the micro ammeter . We have always chosen to 
carry a background rate of ^2 to ^8 counts/second so 
that we were sure that the "zero" was not negative. A 
negative zeroing of the microammeter would have given us 
a false (low) value for the integrated beam current. 

Our next problem is that of maximizing the sig- 
nal to noise ratio. Only in the polarized light inten- 
sity scalar is the S/N ratio low enough to merit concern 
The only option, after external noise sources have been 
eliminated, is to try adjusting the discriminator of the 
single channel analyzer. We have found no appreciable 
changes in S/N over a wide range of discriminator set- 
tings. (The S/N ratio is a function of the pressure and 
ranges from ^3 to ^100.) We have therefore arbitrarily 
chosen a discriminator setting which results in a con- 
venient count rate. Again, see Pegg (1970) for further 
discussion. 

3.4 Taking Data 

The steps followed in the actual data taking 

are j . - 

1. Warm up phototube and electronics for ^12 
hours . 



2 . Warm up accelerator for vl hour 

3. Turn on beam in center port, adjust current 
and energy 

4. / Steer beam into left port, adjust current 

and energy 

5. Go downstairs and . . . 

6. Start flow of He into target chamber, adjust 
pressure 

7 . Assume data taking position at end of left 
port beam tube 

8. Reset all scalars 

,9. Switch beam -off by remote switch 

10. Start scalars for background count 

11. Record values on scalars 

12 . Turn beam back on 

13. Rotate polarizer (by hand) for passing E 
vector vertical 

14. Reset scalars 

15. Start scalars 

16. Record values 

17. Rotate polarizer 90° (by hand) 

18. Start scalars 

19. Record values 

20. Repeat 13-19 with occasional background 
runs (beam off) 

A normal data run will consist of first a back- 
ground reading for the polarization and normalization 



scalars; then 7 sets of "polarization vertical" - 
normalization - "polarization horizontal" - normalization 
then a background run; then 6 sets of data; then a back- 
ground run; then 7 sets of data; and finally a last back- 
ground run. Typically, a data run (20 measurements + 

4 backgrounds) will take ^3/4 hour exclusive of tune up 
time for the Van de Graaff (another ^1/2 hour). Table 4 
is an actual data page from the laboratory notebook. 
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TABLE 4 


DATA FROM THE LAB NOTEBOOK 
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SECTION IV 
ANALYSIS 


4.1 Calculating tt 


The linear polarization fraction tt is defined 


as 


TT* 


*// 


4.1.1 


where I /; and 1^ are the intensities of the light with 
electric field vectors respectively horizontal (parallel) 
and vertical (perpendicular). These intensities are 
determined from the number of counts showing on the 
scalars at the end of each parallel and perpendicular 
data acquisition period. We must of course subtract any 
background counts which are included in the scalar dis- 
play. Furthermore we must normalize these intensities ' 
to integrated beam current (as discussed in section 3) . 
Therefore the intensities are expressed as ... 


X, 


Q - 8 * 


- 8 , 




Cj. - f 3 z 

- g„ 


4.1.2 
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where in both cases the numerator represents signal 
counts, i.e. total counts (on the light intensity scalar) 
minus background counts. The denominator in both cases 
is proportional to the number of protons which have 
passed through the target chamber, represented by total 
counts (or the normalisation-: .scalar^... minus, background 
counts. Thus \f ( and 'Ij_ hire the intensities per proton, 
and are -therefore independent of small fluctuations of 
the proton beam current. We then have the following 
definitions . . . 

c //' c x = The number of counts on the light in- 
tensity scalar for the parallel and perpendicular data 
acquisition periods. 

B x = The average number of background counts 
on the light intensity scalar. 

N // ' = The number of counts on the normaliza- 

tion scalar for the parallel and perpendicular data 
acquisition periods. 

B N ~ The avera 9 e number of background counts 
on the normalization scalar. 

At first glance a non zero is not expected, however we 
purposely set the zero adjustment of the Keithly micro- 
ammeter to give a small background on the normalization 
scalar. This was done so that we could be assured (by 
checking background counts ) that the microammeter zero 
adjustment was not set negative. Therefore we were 
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assured that no counts were lost. 

We then do the actual calculation of the linear 
polarization fraction tt from ... 



- g* ___ C.J. - 8r 

//„ - 

ZJj u - Sr 

N >, “ - 0 V 


4.1.3 


Note that a calculation of tt requires two data acquisition 
periods, one for the parallel and one for the perpendicular 
orientation of the polarization analyzer. A typical data 
run then consisted of twenty measurements of I /y and I L , 
and then the twenty calculations of tt. The average back- 
grounds Bj and were found from four measurements taken 
before, during and after the data run. See again table 4. 

The actual calculations are done on the University 
of New Hampshire’s IBM Call 360 time sharing computer 
system. 


The computer program which does this calculation 
(called NEWPOL and listed in Appendix C.2) will print out 
C m N //f Cj. , Nj_ , and tt all twenty times and will then print 
a mean tt , the standard deviation, and the standard error, 
of the mean, the beam energy, current, the He target pres- 
sure, and the number of tt ' s deviating by more than two 
standard deviations from the mean tt. Table 5 shows, the 
program output created from the data of table 4. 



TABLE 5 


OUTPUT OP THE PROG BA!! 


NEW? CL 


LINEAR POLARIZATION FRACTION 
# OF. RUNS TO ANALYZE IS?1 
LIST ALL P?YES 


V 

C 

H 

C 

P 

899 

1467 

818 

1 502 

-0.06165 

902 

1523 

815 

1533 

-0.05680 

863 

1533 

852 

1 530 

-0.00874 

782 

1486 

838 

1 524 

0.02395 

792 

1531 

795 

1546 

-0.00290 

*809 

1536 

776 

1553 

-0.02561 

809 

1558 

780 

1553 

-0.01775 

84 7 

1430 

795 

1514 

-0.06218 

894 

1559 

784 

1565 

-0.071 26 

804 

1569 

784 

1 546 

-0.00594 

819 

I 549 

802 

1 547 

-0.01047 

824 

1550 

873 

1553 

0.02955 

900 

1564 

793 

1567 

-0.06777 

865 

1633 

883 

1638 

-0.00273 

952 

1634 

838 

1557 

-0.04292 

906 

1 587 

876 

160 0 

-0 .02064 

871 

1 604 

888 

1606 

0.00957 

892 

1 609 

844 

1 607 

-0 .02856 

936 

1 606 

927 

1599 

-0.00288 

661 

142! 

880 

1484 

-0.01030 

300 KEV, 

-20 U 

HE, 

11 .5 U 

A H+, 20 TRIALS 

ST DV « 

-02991 

0 PTS DV > 

2 SIC S.E.= .00669 

<P-SE) =- 

• 02849 

MEAN 

P = — 0 i 

• 02 1 80 CP + SE>=-.02 5! 1 
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4.2 Polarization Measurement Error Analysis 

There are of course two types of error inherent 
in this experiment. Random errors occur due to the nor- 
mal variations (fluctuations) in counting rates? and 
systematic errors occur if faulty equipment and/or 
methodology. is used. The random errors present a problem 
only in that they limit the precision of the calculation 
of 7 t; any systematic error could reduce the accuracy of 
the experimental result. 

We will first discuss the random errors present 
in this work. The linear polarization fraction tt is 
calculated from Eq. 4.1.3. The variables which. appear in 
this equation are simply the numbers which are displayed 
by the scalars at the conclusion of each of two 1 1/3 
second intensity data acquisition periods, as well as 
the average of four 1 1/3 seconds background count 
acquisition periods. Every one -of these scalar displayed 
numbers is subject to normal statistical fluctuations, 
thus when tt is calculated, the precision of a single cal- 
culation is quite poor. We therefore repeat the measure- 
ments from which tt is calculated twenty times for each 
"run," and . take many runs for every "final" value of tt. 

The 7T. which we will report as a result is the average 
of Nx20 values, (where N is the number of runs at some 
particular pressure, beam current, and beam energy). Thus, 
the precision of this work is enhanced by repeating our 
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measurement a multitude of times. 

- order, to. analyze our random errors, we must 
take care to distinguish between sample variables and 
population variables. We will follow the notation of 
Parratt (1971) and Wilson (1952) where the sample (run) 
mean , standard deviation, and standard error are sym- 
bolized by m,s,.and s m ; while the population mean, stan- 
dard deviation, and standard error are symbolized by y, 
a, and a m . We are of course most interested in y and 
a m s: *- nce We will report the linear polarization fraction 
in terms of our best estimate of the parent population 
mean, and its standard error* The procedure is outlined 
in the following formulae where tt^ is the calculated tt 
from one set of measurements for the j th run (i.e. 1/20 
of a run) and n=20, All tk ^ are equally weighed. For 
each run, where j is the run index 


"V ■ - %■ 


4.2.1 


while for the best estimate of the parent population 
mean ... 


N n 

K-j 


4.2.2 


where N is the number of runs at some particular He target 
gas pressure, beam current and energy. It now remains to 
calculate and estimate the dispersions in the sample and 
parent populations. The standard deviations s and a are 
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suitable dispersion indicies. For each run we have 



4,2.3 


while for the parent population we have 



4.2.4 


However, p is never known exactly, only estimated as u f 
from Eq. 4,2.2. The best estimate of a {see the appendix 
of Bacon (1953) is then 



4.2.5 


Since in this work we take many runs and calculate 
a mean n for each, it is useful to also calculate the stan- 
dard error a m (also called the standard deviation in the 
mean). The experimental standard error S m is 


s„ = 


1 

' E. , / N . vl 1 

1 T («y -%-a «y) J 


4.2.6 


where - Wj 


is the grand mean, or mean of means. 


s m is more easily estimated from the standard deviation s 

by 




n> 


s 


4.2.7 



where in our case n=20, the number of r> ’ s in one run. . We 
are in fact more interested in the standard error cf 
the parent distribution. 

4.2.8 

and 




L isr> 

A I'poc 


<r. 


07 


jTTv 7 


cr 


4.2.9 


Again u is unknowable exactly and we must estimate 
Using 4.2.9 with 4.2.5, we arrive at the useful result 
that 



4.2.10 


According to statistical theory, the mean tt of 
any run has a probability of 68.3% of falling between the 
values of y ± 

We calculate and report as follows . . . for 
every run of n=20 tt's, we calculate the mean for the run 
m_. {Eq. 4.2.1); an estimate of the parent population 
standard deviation for the run ct^ (Eq. 4.2.5, with N = 1) 
and the estimated population standard error a m (Eq. 4.2.10, 
with N = 1) . 


When sufficient runs at a particular He target 
gas pressure, beam current and energy have been' 



accumulated, the Nx2Q tt values are combined to estimate 

the population mean...p\ {Eq. 4 .2. 2) usually written 

simply as ir ; the population standard deviation a 

{Eq. 4.2.5); and the population standard error a 
■ m 

(Eq. 4.2.9). 

In the first case where single runs are analyzed, 

the computer program NEWPOL is used (Appendix C.2), for 

calculation of the accumulated result the computer program 

STAT is used (Appendix C.3). 

When our results are given in the next section, 

we will report the estimated population mean and the 

standard error a . 

m 

We must now discuss possible systematic errors 
since this type of error will reduce the accuracy of any 
experiment. Some typical examples of systematic error 
are . , . 

1. Subconscious bias on the part of the 
observer 

2. Incorrectly calibrated instrument (s) 

3. Misaligned apparatus 

4. Lack of correction for changing 
environmental conditions 

False assumptions of independence of result 
from an experimental parameter 
6 . etc . 

In this work, the greatest possibility for a 
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systematic error lies in whether or not our polarization 
detector is biased. We evolved a simple method for 
determining the existence of any bias. Our polarization 


detector . housing was modified such that it was possible 
to rotate the entire detector, about its line of sight, 
by 90° . When we then compared it results arrived at for 
the detector normal (N) position and the results for the 
detector rotated (R) position, we could unfold any de- 
tector bias. We chose to do this procedure over a range 
of proton beam energies with constant beam current and 
He target gas pressure. Many runs were taken at each 
energy with detector (N) and (R) . The means tt xt and tt^ 
were calculated at each energy along with the estimated 
population standard deviations c x1 and c n . We are then 
able to use the "student" t test for the significance of 
the difference of means as described by Spiegel (1961) . 
The t test is derived from small sampling theory 
(although, under conditions of large samples, it becomes 
equivalent to the z statistic) . Using our notation . . . 


t -- 


sty tj ~ * 


4.2.11 


n % 


/> At, 


where n again is 20 and (N R ) is the number of runs 
with detector N(R) . Also, 


CT r 


A* + n Sp 


n 


n N, 


4.2.12 


N 


- 2 
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where S N (S^) are the detector N (R) standard deviations 
calculated from 


5 - 


‘ 4 / Al 
a:j a-.i 


N 



4 . 2.13 


Note that 



where a was defined in Eq. 4.2.5. 


Having done these calculations, we are then faced 
with the statistical decision making process. In order to 
compare two means, we must do a two-tailed test on the t 
statistic. Let us state the hypothesis H Q : there is no 

significant dif ference • in the means? and : there is a 

significant difference in the means. It is customary to 
make the test at both the .05 and .01 levels of signifi- 
cance. The range of acceptable t values at the ,01 and 
.05 levels is shown in table 6 for various degrees of 
freedom *7 = n +nN R -2 .If when we calculate t 
from Eq. 4.2.11, t lies outside the range given by 
table 6 {for a particular *7 and level of significance), 
we would have to reject H 0 at* that level of significance.. 

A rejection at the .05 level means we have a 5 percent 
probability of having made a "false" rejection, and a 
rejection of H Q at the .01 level implies a probability of 
1 percent of a "false" rejection. If we have a rejection 
at .05, but not at .01 then E 1 is "probably" true, i.e. 
there may be a significant difference in the means. 

A computer program TTEST (Appendix C.4) was used 
to calculate the t values from Eq. 4.2.11 for both N and 
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TABLE 6 

T-TEST SIGNIFICANCE LEVELS 
SIGNIFICANCE LEVEL 



.01 

.05 

20 

-2. .84 < t < 2.84 

i 

-2.09 < t < 2.09 

40 

-2.70 < t < 2.70 

-2.02 < t < 2.02 

60 

-2.66 < t < 2.66 

-2.0 < t < 2.0 

120 

-2.62 < t < 2.62 

-1.98 < t <1.98 

00 

i 

-2.58 < t < 2.58 

-1.96 < t < 1.96 


R orientations of the detector. The results are shown 
table 7. Only at 450 keV proton beam energy is there 
any possible significant difference in the means . The 
last entry in table 7 compares the mean of the (N+R) 
results with the mean of the N results, and shows no 
significant difference. We therefore conclude, that our 
polarization detector has no inherent polarization form 
bias and that said detector does not contribute any 
systematic error to our experiment. We will henceforth 
combine the results of N and R runs taken at the same 
pressure, energy and current. 

A second possible source of systematic error in 
this work could be caused by any drift in our electronic 
which was constant in sign for the greater part of a run 
In order to average out any such drifts, the order of 
data taking was varied during a run. The linear polar- 
ization fraction was always calculated for time adjacent 
I// an< 3 ij. values. Our data taking order would be . . . 

[(li.Ni) (I,, HP ] [ (I// ) (Ij.Nl ) ] td^N, ) (I # N,,) ] ' 

. . . etc. where the rectangular brackets denote a cal- 
culated 7r. It is our contention that any "constant" 
drifts would be averaged out by this method of data 
taking. Any random drifts would of course average them- 
selves out of our calculations. 

Another possible source of systematic error lies 
in the fact that stray magnetic fields can cause preces- 
sion of the atomic electrons and, therefore, depolarize 
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TABLE 7 

TEST FOR SIGNIFICANCE OF DIFFERENCE OF MEAN tt 
FROM NORMAL AND ROTATED DETECTOR 


BEAM 

ENERGY 

t VALUE 
N Vs R 

# of DEGREES 
Of FREEDOM 

DIFFERENCE IN MEANS 
SIGNIFICANT 




.15 

,01 

(keV) 





150 

0.00 

138 

No 

No 

200 

1.75 

118 

it 

it 

250 

* .25 

118 

H 

' 

ii 

300 

-.48 

238 

M 

» 

350 

.96 

118 

II 

ii 

400 

1.76 

118 

II 

ii 

450 

2.48 

138 

YES 

it 



t VALUE 
(N+R) Vs N 


i 


450 

-1.23 

218 

No .*■ 

No 
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the light from a collection of atoms. Feof i.lov (1961) 
discusses this problem and gives the results of Breit 
as follows ... 



and 4.2.14 

where ir 0 is the zero magnetic field linear polarization 

fraction, it is the polarization fraction measured at 

magnetic field strength H, t is the mean life of the 

transition, g is the Lande g factor and <p is the angle of 

rotation of the plane of polarization. Using . . . 

e = 4.8xlO~ 10 stat c 
m = 9.1xl0~ 2 0 gm 
c * ' 3x10 1 0 cm/sec 

estimating t*2x 10” 9 seconds (from Dotchin et al, ) and cal- 
culating g for the 3*P states from 

, 4.2.15 

where j=l, s=0, and i- 1 so that g=l. We then have . . . 

7 Z = 7rJJ tft.ois)"]' 1 4.2.16 

Measurements of the magnetic field strength H in 
the vicinity of the target chamber were made using a 
HeliFlux Magnetic Aspect Sensor type RAM-3. These measure- 
ments were made with the Van De Graaff steering magnet 
energized. In no case were we able to find a field in 
excess of ^3/4 gauss. Using this in Equation 4.2.16, we 
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get . . . 


771 - V ( i +. 0007) 


and conclude that the magnetic depolarization is negli- 
gible when compared with the random error present. Simi- 
larly we conclude that $=*0 and the plane of polarization 
is not rotated by the magnetic field* 

We must also concern ourselves with collisional 

depolarization. Collisions between helium atoms in the 

target chamber, as well as with the chamber wall, will 

tend to equalize the magnetic substate populations and, 

therefore, reduce the polarization. Referring to a 

mean free path calculation done in the next section 

(4.3), we find the mean free path (using the geometric 

cross-section) of the helium atoms to be =* ,47cm. The 

rms velocity of the helium atoms is . '\T me = f 3 K ~r ] 

v ~~m~ 

which for room temperature helium works out to 

s ~ /. y x /o* c/s/ sec • Therefore the mean collision 

time is 3.3x10 4 sec. However, we must also account 
for the chamber walls where the radius is ^2. 5cm. Thus, 
we can expect a wall collision in t^l.8xl0~ 5 seconds. 
Since even this time is 4 orders of magnitude larger 
than the mean life of the 3 *P state (1.8xl0“ 9 sec), we 
can ignore the possibility of depolarizing collisions. 

Another mechanism which could cause a systematic 
error is that of cascading, where the 3*P level is popu- 
lated not only by direct proton excitation, but also by 
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downward transitions from n'S and n'D levels (n>3) . In 
general, cascading into the 3*P level from a number of 
different upper levels will tend to equalize the mag- 
netic substate populations and, therefore, reduce the 

amount of polarization exhibited by the 3 1 P * 2 l s 

transition. Heddle (1962) has shown that if the per- 
centage of the 3*P population due to cascading is C% 
then the experimentally observed polarization must be 
multiplied by the factor (1 + c/100) in order to correct 
for the cascade depolarization. C is most easily cal- 
culated following the method in appendix b of Van Eck f s 
(1964) paper, using the transition probabilities listed 
i n G a ktiel (1960). However, in order to do this calcula- 
tion, we need the values of various experimental emission 
and excitation cross sections. All the necessary experi- 
mental cross sections for this calculation are not yet 
in the literature. Since the agreement between theoreti- 
cal and experimental cross sections is so poor, we have 
decided not to do the calculation with the theoretical 
values, but to instead rely upon the experimental deter- 
minations of cascade given by Van den Bos (1968) and by 
Thomas (1967). In both of these works, the authors use 
their own unpublished data to calculate the cascade per- 
centage c, which they then report. Van den Bos reports 
values for c of 8% at 100 keV and 6% at 150 keV? while 
Thomas reports 4% over the energy range 150 - 1000 keV. 

We will, therefore, correct for cascade as follows: 
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at 100 keV use 8%, at 150 keV use 5%, and above 150 keV 
use 4% for the cascade correction. 

We must next discuss the effect of resonance 
trapping on the polarization of the emitted light. 
Resonance trapping occurs when the beam excited 3 : P 
level decays to the 1 1 S level and the resultant photon 
is trapped by a ground state He atom. The ground state 
atom is then of course excited to the 3 l P level and may 
then undergo a 3*P — -*■ 2 J S transition. This event re- 
sults in our detector seeing light from atoms which were 
not space quantized by the proton beam, and, therefore, 
results in an experimental polarization which is lower 
than the free atom value. In order to estimate the 
amount of trapping which occurs, we turn first to the 
work of Lees (1932) and Lees et al. (1932) . These works 
used photographic techniques to determine the excitation 
cross sections for various He levels excited by electron 
impact. When Lees' photographs showed X 5016 A light 
emanating from He gas which was not impacted upon by the 
beam, he and Skinner concluded that 3 1 P-1 1 S radiation 
from the beam area was being trapped and a portion of it 
would then be emitted as 3 1 P-2 1 S light. They showed that 
the proportion of light absorbed = f - ^ T< where n is 

the number density, r is the target chamber radius, and 
T a ^he atomic absorbtion coefficient. 

Their values for were 2.3xl0“ 15 (experimental) 
and ^92x10 1 (theoretical). Once we know the 



percentage trapped in the 3 1 ? level, the ratio of the 
transition probabilities A (3 l P-2 *S) /A ( 3 1 P-1 1 S ) - .023 
will give us some indication of the percentage of 
3 1 P-‘1. 1 S which is first trapped and then radiated at 
3 1 P~2 1 S radiation. Of course, the reradiated 3 1 P-1 1 S 
radiation may be retrapped, etc. Table 8 shows the re- 
sults of these estimations as a function of pressure. 
This issue is obviously not resolvable with this data,. 

Recently more comprehensive work on the imprison 
ment of Resonance Radiation in gases has been done, 
notably by Holstein (1951) and Phelps (1958). Holstein' 
mathematics are not suited to our low pressure work, and 
in fact are not useable, Phelps on the other hand using 
Holstein's concept presents a transmission probability 
graph from which the transmission of the radia- 

tion, of He (in terms of target chamber radius) can be 
abstracted. This transmission probability is then con- 
verted into % trapped (see table 8). We note that the 
Phelps results fall between the two extremes of Lees’ 
work. Using Phelps' results in the worst possible case, 
where 93% of the 3 1 V-l l S radiation is trapped, and 
applying the factor of .023 from the ratio of transition 
probabilities, we find a trapping correction T of 2.1%. 
This has assumed single trapping events. In order to 
account for multiple trapping we arbitrarily raise that 
figure by ( . 023) x ( • 5) - . 012 and estimate the trapping 
correction T as 3.2%. Estimated trapping corrections 
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TABLE 8 


TRAPPING CORRECTIONS 


He 

PRESSURE 

mtorr 

LEE’S ' 
EXPERIMENTAL 

% TRAPPED 

LEE ’ S 

t THEORETICAL T 

a a 

PHELPS 

RESULTS 

TRAPPING 
CORRECTION 
T in % 

.01 

.15 

5.9 

5 

.17 

.05 

.75 

26. 

20 

.69 

.2 

‘3.0 

70 

45 

1.6 

.6 

8.7 

97 

80 

2.8. 

1.0 

14 

99 

90 

3.1 

1.5 

20 

99 

93 

3.2 
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are listed in table 8 for the various He pressures we 
have used. We can then correct for resonance trapping 
depolarization by multiplying our experimental tt by 
(1 + T/100) . We will later discuss another type of 
trapping which may cause error, namely trapping at the 
2 1 S level of a 3 1 P-2 1 S photon. 

* 

Since we are already making a cascade correction 
of the form (1 + C/100), and rather than compound the cor- 
rection, we define D = C + T as the unpolarized light 
contribution to the 3 1 P-2 1 S radiation. Our correction 
therefore takes the form tt = tt (1 + D/100) . Table 9 

G Xp 

shows D tabulated for various proton beam energies and 
Helium gas pressures. Generally this correction is 
smaller than the standard error of the mean of tt. 

We must also discuss the effect of the He gas 
purity on the measured polarization. Since we work at 
quite low pressures, 1.5 mtorr and below, the gas density 
is quite low and as stated earlier the mean collision 
time is roughly 3x10” 4 sec. However the 3 3 P state has 
a mean life of 1.8x10 9 sec (Dotchin — private communica- 
tion) and decays of that state should be extremely free 
of effects from the very infrequent collisions with not 
only other He atoms, but also with the rare impurity 
atom. We have used two sources of He in taking data, the 
first was standard grade ^98% pure, and the second was 
research grade 99.999% pure. The results obtained were 
identical (within S.E.) and clearly impurities pose no 
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TABLE 9 

ESTIMATED UNPOLARIZED LIGHT CORRECTIONS 





He 

Pressure 

(m torr) 





.01 

.05 

.2 

.6 

1.0 

1.5 


100 



9.6 





150 

5.2 

5.7 

6.6 

7.8 

8.1 

8.2 


200 



5.6 




Proton 

Energy 

250 



5.6 




(keV) 

300 

4.2 

4.7 

5.6 

6.8 

7.1 

7.2 


350 



5.6 





400 



5.6 





450 

4.2 

4.7 

5.6 

6.8 

7.1 

7.2 
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problem to this study. 

A simple problem which could have caused a gross 
systematic error was dust. Any substance which might 
adhere to our ...UN- 32 polaroid analyzer and scatter light, 
could cause an artificial difference in the intensities 
I l{ and as the analyzer was rotated between it's two 
positions. We therefore took great care to keep the 
polarization analyzer clean and dust free. To insure 
that an error caused by dust would be reduced, we spread 
our data taking over many days and many analyzer clean- 
ings, assuming that in the long run unescapable small 
errors would average out. 

Another possible source of error is due to the 
shape of our proton beam. It is roughly rectangular in 
cross-section, ^2x10 mm. This is quite a large beam area 
for an experiment of this type, and we need to be con- 
cerned about those processes which might occur within the 
beam region. We are forced to accept this unusual geom- 
etry due to both the idiosyncrasies of our electrostatic 
focusing system and the experimental needs of Dotchin 

sli.* The beam region in the target chamber is a plasma 
region, and a Faraday rotation of the plane of polariza- 
tion of light is conceivable. Faraday rotation occurs 
due to birefringence in the plasma, the birefringence in 
turn is due to a magnetic field (e.g. the earth's field). 
The birefringence also means that the indices of refrac- 
tion are different for right and left circularly 
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polarized light. Since a linearly polarized wave can be 
considered to be the superposition of two counterrotating 
circularly polarized waves, and since the indices of re- 
fraction are different for each component (hence the phase 
velocities are different) , the plane of polarization of 
the linearly polarized wave is rotated. The amount of 
rotation depends upon the differences in the two indices 
of refraction, and the thickness of plasma penetrated by 
the wave. An outline of the Faraday effect is found in 
a text by Marion (1965) and a very good generalized 
treatment is available in the work of Larson (1967). The 
two indices of refraction are found from 

w(w*s) ] 4.2.1 

T* V )/ A/e "7 

where o 0 is the plasma frequency uV - J rn ' N e 

is the electron density; s is the gyrof requency of the 

* ' ' ' - D 

electrons in the magnetic field $ = * B is 

the magnetic field component along the propagation direc- 
tion; and oj is the frequency of the light from the 
3 1 P~2 I S transition. The rotation of the plane of polari- 
zation is given by 


ciJL 
d L 



4.2.18 
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' iSl . 

WnerS JT 1S the rotat ^ on i n radians/cm of penetration. 

The important quantity is obviously (n + -n_) this' in turn 

depends upon w±s. For our work we have - Z. 76 *'<>'* 

o 

and assuming B = . 5 gauss, s == 8.78xl0 6 . Using formula 
4.2.17 we find that the difference between n + and n_ is 
on the order of one part in 10 15 or (n + -n_) * 10“ 15 and 

therefore =* 6x10"- 11 radians/cm. We can therefore 

state that any Faraday rotation of the plane of polariza- 
tion of the light from the He atoms is totally and com- 
pletely ignorable. 

Another problem which could have caused a sys- 
tematic error was caused by soft x-rays from the Van de 
Graaff tank. These x-rays contributed to the background 
count rate for our phototube. However to take a back- 
ground count measurement, we shut off the belt charge 
and therefore the beam and x-rays. Thus we had been get- 
ting false (low) background count measurements. This 
problem has been solved by shielding our phototube with 
lead bricks so that essentially no x-ray counts are de- 
tected. 

We have taken pains to eliminate other sources 
of systematic, error. The detector has always been care- 
fully aligned at 90° to the beam direction with the use 
of a "square." We have checked our electronic 1 s count- 
ing rates by feeding a simultaneous pulsed signal to all 
systems and observing identical numbers on all scalars. 
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We, therefore, conclude that systematic errors 
in the polarization measurement have been reduced to 
the extent that they are not significant in comparison 
to the random errors present. 

4.3 Other Parameters 

Since it is our aim to find the linear polariza- 
tion fraction as a function of Helium target gas pres- 
sure, proton beam energy, and proton beam current, we 
must both measure (and estimate the error in our measure- 
ments of) these parameters. 

We measure the Helium target gas pressure with an 
instrument which, according to the manufacturers specifi- 
cations is capable of five digit accuracy. This MKS 
Barratron capacitance manometer will not in of itself 
be a source of error. We could, however, have a pres- 
sure error if there is "significant" outgassing in the 
piping which connects the target chamber to the pressure 
head. See figure 7, This piping is primarily 1/4" O.D. 
copper tube soft soldered to brass and/or stainless 
steel fittings. In our observations of the behavior of 
this section of the vacuum system, we are able to state 
that under various valve closed conditions, the pressure 
indicators available have not shown excessive outgassing. 

Another possible source of error in the pressure 
measurement system is the fact that the pressure measure- 
ment port of the target chamber is located %2.5'cm 
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"upstream" of the observation window, see figure 6. We 
must then ask whether or not the He gas is undergoing 
Laminar flow and exhibiting a relatively large pressure 
gradient in the target chamber, or whether the He is 
streaming through the chamber with little pressure 
gradient. The latter case is of course the desirable 
situation. We can decide between cases by doing a simple 
mean free path calculation. We have for the mean free 
path l, Young (1964) . . . 



K T 

y/T ir r 2 P 


4.3.1 


where k is the Boltzman constant k=-l . 38x10" 2 3 j/°k , T is 
the temperature, estimated 300°k; rr 2 is the geometrical 
cross section estimated as 3xlO” 20 square meters; and P 
is the pressure in N/m 2 , We then have . . . 



• 0CfV j 


4.3.2 


The largest pressure at which we work is 1.5x10“ 3 torr = 
1.5x10 3 mm of H^=,2 N/m 2 . Therefore, the mean free path 
for our highest pressure is £=.47m=47cm. We then con- 
clude that the Helium gas is streaming through our cham- 
ber and that no significant pressure gradient exists. 
Therefore, the fact that our pressure measuring port is 
upstream of the observation window will cause no signifi- 
cant pressure error. 
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Our next concern is whether or not the proton 
beam energy at the target chamber is the same as at the 
accelerator. We have previously discussed, in 
section 2,2, the accuracy of the beam energy readout of 
the Van de Graaff accelerator system. Our estimate of 
the accuracy of that readout was ^8 keV. We must also, 
however, ask whether or not the proton beam energy is 
effected by the presence of the target gas, or "back- 
ground" gases in the beam tube. We can use the 
stopping power data compiled by Allison (1953) to cal- 
culate dE/dZ, i.e. the proton beam energy loss/cm. In 
the energy range of 100 to 500 keV, the greatest stopping 
power of Helium for protons (occurring at 100 keV) is 
7.3x10" 15 evxcm 2 /atom. We can easily calculate the num- 
ber of atoms/cm 3 from the ideal gas law . . . 


PV = nKT 


4,3.3 


where P=1 . 5x10" 3 mm of H g =22 dynes/cm 2 , K is the Boltzman 
constant 1.38x10 16 erg/°K; and T is estimated room 
temperature 300°K. Therefore, the number density of He 
atoms in our target chamber is £ = ± s : y/ ? * /0 ' J 

Y n ' 


Multiplying the stopping power by the number density, we 
arrive at the energy loss/cm = 35xlO“* 2 ev/cm. Thus the 
beam energy lost to the target gas is completely negli- 
gible. Since the "background" gas pressure is ‘MO 2 less 
than the target pressure, we can estimate an Energy loss 
to "background" gas of .0035 ev/cm, again totally 



negligible. We, therefore, conclude that the systematic 
error in the beam energy is not significant when com- 
pared to the v8 keV instrumental uncertainty. 

The next important parameter is the proton beam 
current, which is detected within our Faraday cup (dis- 
cussed in section 2.3). We have avoided the common 
secondary electron loss problem in the cup in two ways; 
the proton beam impacts on Vycor glass rather than on 
metal thus producing fewer secondaries; and we have 
lengthened our cup and provided a small entrance aper- 
ture such that secondary electrons could only escape 
through a very small solid angle. Considering other 
factors such as leakage through insulators, etc., we 
estimate an overall accuracy of v5 percent on the collec- 
tion of beam current by our Faraday cup. 

The readout of the proton beam current measure- 
ment is then taken two ways for two purposes. The 
Keithly microammeter meter readout is used to tune the 
Van de Graaff and to provide a number for the comparison 
of the linear polarization fraction with beam current. 

We estimate an instrumental accuracy of v2 percent. 
However, in practice, the beam current is a very un- 
stable quantity, if for example the Van de Graaff has 
been tuned to give a vlo U A H + beam in the left port, the 
microammeter will show short term fluctuations in the 
worst case of from v 9 to *llpA. The reported beam cur- 
rent of viOyA is the observers best estimate of the time 



84 


averaged beam current, uncertainty v5 percent. This 
problem is further complicated by drifts in the Van de 
Graaff’s output current. During the course of a run 
(^45 minutes) the beam current may change, e.g. from vio 
to ^12uA. We will, therefore, state an overall average 
beam current uncertainty, for any run of ^<15 percent. 

The second type of beam current readout is taken 
from one of our scalars and used to normalize light in- 
tensity to beam current. This scalar readout is in 
arbitrary units since an absolute number is not required, 
see Eq. 4.1.3. Our main concern here is repeatability, 
i.e. a given number of protons impacting on the Vycor 
end window of the Faraday cup should produce the same 
number on the normalization scalar. We have no direct 
way of confirming this repeatability, however, our ob- 
servations of the performance of this equipment has en- 
abled us to conclude that so long as the Keithly micro- 
ammeter and other electronics have been warmed up for 
hours prior to a run, the repeatability will be quite 
good. Once the electronics have warmed up, the random 
errors on the normalization scalar far outweigh any re- 
maining systematic error. 

4.4 Data Manipulation 

From the time data is first collected and written 
in our lab notebook, until the time a. "final" 'ir is re- 
ported, much, computer manipulation occurs. We must 



85 


analyze each and every run prior to joining the results 
of that run to previous results (to guard against major 
systematic system error) . We must maintain a separate 
file for every target gas pressure, proton beam Energy, 
proton beam current, and detector orientation (to en- 
able us to separately analyze it as a function of each of 
those parameters). Table 10 shows the number of (20 
polarization measurements) runs taken at each beam en- 
ergy and target gas pressure for both N and R detector 
orientations and for normal and low beam currents, a 
total of 116 runs. Due to the large amounts of data in- 
volved, we have established a system whereby all of the 
data and calculated tt's are stored in the University of 
New Hampshire's computer (call 360 system). 

All the data taken for one run (Table 4) will 
be entered into the computer memory, as a program, using 
the data statement. Table 11 shows the section of a 
data program relevant to the data of Table 4. All of our 
data is hand typed into the computer memory and onto 
punched paper tape. We have used a consistent naming 
system for these data programs. They are labeled as 
DATAXXYY, where XX is a pair of numbers representing the 
month, and YY is a pair of numbers representing the day 
of the month. Any data program will include all the 
runs done during one day. These data programs contain 
nothing but data statements and are useless unless merged 
with an operational program. 
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TA3LE 11 


DATA P ROC- PA- ' 



datao4 i 4 


9:36 AUGUST 27, 1973 


REM DATA PI 
REM ENERG 
DATA 300, 
REM DATA. 
DATA 899, 
DATA 7 62, 
DATA R09, 
DATA 604, 
DATA 900, 
DATA 906 , 
DATA 936, 


4 1 4 

V# PRESSURE, CURRENT 
•2,1 1 • 5, A, 20, 49, 8, 
..COUNTS VERT, NORM 
1467,818,1502,902, 
1486,838, 1 524,792, 
1558,780, 1 553,847, 

1 569,784, 1 546,81 9, 

1564,793,1567,885, 
1587,878, 1 600,871 , 
1606,927,1 599,861 , 


,4, 20, BACKGROUNDS 

35,8,52,9,47,10 

,HOR I 7, NORM, VERT, NORM. . .ETC 

1523,815,1533,863,1533,852,1539 

1531 .795. 1546.809. 1538.778. 1 553 

1430.795. 1 514,8 94, 1559,7 84, 1 565 

1 549.802. 1 547.824. 1550.873. 1 553 
1 633, 883, 1 638,952, 1 634, 838, 1 

1604.888. 1606.892. 1609.844 . 1 
1421,860,1484 



o t.n 
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The first operational program used in NEWPOL 
(appendix C.2). Its function is, when merged with 
DATAXXYY; to separately and sequentially, analyze each 
run; to calculate all 20 irr 1 s for each run; to calculate 
a mean, it and the standard deviation and error; and to 
print for each run the raw data and the calculated re- 
sults. See again Table 5. 

In order to combine runs, their individual results 
have to be cumulatively stored. We have a series of pro- 
grams which will analyze the runs within DATAXXYY and, 
instead of printing the results, add to storage files the 
20 tt * s from each run according to target gas . pressure , 
beam energy, and detector normalization. These programs 
are labeled PUTPOL, PUTPOLR, and PUTPOLP (appendicies 
C. 5,6 and 7). The files into which the results are placed 
were created from the keyboard by user commands and ini- 
tialized by using PUT (appendix C.8). PUTPOL is merged 
with detector N data with Normal beam currents. It will 
only analyze and file .2m torr runs at energies from 
100 to .450 keV (by 50 keV steps). These files are la- 
beled PXXX where XXX is the beam energy in keV. PUTPOLR 
resembles PUTPOL except that it is used on detector R 
runs and the files are labeled PXXXR. Finally PUTPOLP 
is used for. all beam energies and target gas pressures 
excepting ,2m torr . This program was put into use 
after we had concluded that our experimental result was 
independent of detector orientation, therefore, the 



output files were labeled only by beam energy and target 
gas pressure. We selectively analyze and file only 
normal beam current data with PUTPOLP. The output files 
are labeled T150PYY , T300PYY, and T45QPYY; or TXXX where 
YY represents the Helium target gas pressure: 

01-.01 mmtorr, 05=,05 mmtorr, 2=. 2 mtorr, 6= .6 mtorr, 
1=1 mtorr, and 15 = 1.5 mtorr; and XXX represents beam 
energy for. 100, 200, 250, 350, and 400 keV where only .2 
mtorr data was taken. See again, Table 10. Low beam 

f 

current (vlpA H + ) data were analyzed and filed using 
PUTPTEST (appendix C.9). The files were labeled T 3 0 OLY Y , 
where YY is defined as above. 

These files are all structured alike as a one 
dimensional array of the form m,l,fr 2 ,ir 2 ,tt 3 , . . . where m « 
the total number of it's (nxN). This form was chosen so 
that we might use IBM's statistical analysis program 
ST^TPACK. Generally, we used our own programs since our 
files contained too many entries for use by STATPACK. 

We, therefore, have over 40 different files rep- 
resenting different conditions of detector orientation, 
beam current, beam energy, and Helium target gas pres- 
sure. In order to analyze the cumulative contents of a 
file, we use the program STAT (appendix C.3). This pro- 
gram will find a mean tr for a file as well as a standard 
deviation, standard error of the mean, etc., see Table 12 
for Typical Output. 

Thus, we have the capability of individually 
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TABLE 12 


OUTPUT OF THE PROGRAM STAT 



STaT 9 S 410 AUGUST 27, 1973 

FILE TO BE ANALYZED IS7T300P05 
U0 P'S IN CALCULATION 
ST DV = #05651 8 PTS DV > 2 SIG 

< P - S E ) = -0.01297 MEAN P = -0.00819 


S.E.a .00478 

(P+SE ) - -0.00342 



analyzing runs using NEWPOL, or of cumulatively filing 
data by parameter for later analysis by STAT. 
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SECTION V 

RESULTS AND CONCLUSIONS 
5.1 Summary of Results 

In this work we have investigated the linear 
polarization fraction of light (X 5016 A) from the 
3*P-2*S transition in Helium as a function of proton beam 
energy, He target gas pressure and proton beam current. 
Because of the large random errors inherent in our experi- 
ment, determinations of tt were made in groups of twenty, 
i.e. twenty measurements for each run. Each run of 
course was done at some particular He pressure, beam cur- 
rent and energy. Any number of runs would be taken on a 
particular day depending upon the behavior of the equip- 
ment and apparatus. The run by run results are listed in 
table 13 in chronological order. Table 14 lists the run 
by run results chronologically within each section of 
beam energy and target pressure. In both of these tables 
we have reported tt in percent, rather than as a (decimal) 
fraction, for the sake of easy readability. 

Table 15 and figure 16 show our results for it 
(in %) as a function of proton beam energy and for con- 
stant He pressure and reasonably constant beam current. 

The graph (fig. 16) uses the corrected it’s. Tables 16, 



RESULTS FOR EACH 


H + 

HE 

ENERGY 

PRESSURE 

KEV 

. MTORR 


100 

.20 

1 50 

.20 

200 

• 20 

250 

.20 

300 

.20 

350 

• 20 

400 

.20 

450 

• 20 


450 

• 20 

400 

• 20 

350 

.20 

300 

.20 

250 

.20 

200 

.20 

150 

.20 

100 

.20 


200 

• 20 

250 

.20 

300 

. .20 

350 

• 20 

400 

.20 

450 

.20 


4 50 

.20 

400 

.20 

350 

.20 

300 

• 20 

300 

• 20 

250 

♦ 20 

200 

• 20 

1 50 

.20 

300 

.20 

300 

.20 

300 

• 05 

300 

• 20 

300 

• 60 

300 

1 .00 

300 

1.50 

300 

2.00 


ABLE 13 

RUN" , CH RONOLOO I C ALLY 


H+ 

C UR.R ENT 

IN 

MICRO-A 

% 

DATA 1216 

4.5 

2.760 

8.0 

1.921 

' 8.0 

1 .507 

8.5 

-0.308 

8.5 

-1 .359 

9.5 

-2.358 

9 . 5‘ 

-3.418 

9.5 

-3.815 

DATA! 2 I 7 

9.5 

-3.477 

9.5 

N -2.194 

9.5 

-3.142 

8.0 

-1 .279 

6.5 

-0.883 

7.5 

0.675 

4.0 

2.624 

5.0 

2.538 

DATA 1224 

10.0 

1 .097 

10.0 

-0.785 

12.0 

-2. 451 

10.5 

-4.498 

12.0 

-4.513 

11.0 

-6.421 

DATA 122 5 

I 1 .0 

-5.832 

10.0 

r4 .585 

1 1 .0 

-3.304 

11.0 

-3.909 

1 1 .0 

-2.224 

10.5 

, -0.991 

1 1 .0 

0.539 

8.5 

1.791 

19.0 

-3.221 

4.5 

-3.358 

10.5 

-1.321 

10.5 

-3.850 

11.0 

-4.766 

11.0 

-4.049 

1 1 .0 

-2.896 

1 1 .0 

-2.196 


S.E. 

*F 

MEAN 


0.569 
0.443 
0.495 
0 .580 
0.730 
0.407 
0.672 
0.585 


0.740 

0.737 

0.640 

0.570 

0.795 

0*7 66 
1 .205 
0.900 


0.505 
0.455 
0.527 
0.669 
0. 529 
0.696 


0.886 

0*604 

0.569 

0.723 

0 .451 
0.430 
0.374 
0.536 
0.300 
0.699 

1 .240 
0.542 
0.175 
0.302 
0.157 
0.107 



TABLE 13 CONTINUED 


H+ 

HE 

ENERGY 

PRESSUR 

KEV 

MTORP. 


300 

• 20 

450 

• 20 

150 

• 20 


150 

• 20 

200 

.20 

250 

.20 

300 

• 20 

350 

• 20 

400 

.20 

450 

.20 

450 

.05 

450 

• 60 


150 

.60 

150 

.05 

150 

.20 

200 

.20 

250 

.20 

300 

• 20 

350 

• 20 

400 

.20 

450 

V20 


300 

1 .00 

300 

• 60 

300 

.05 

300 

.01 

150 

.01 

150 

.05 

1 50 

1 .00 

150 

1 .50 

150 

.01 


150 

1 .50 

150 

1.00 

150 

• 60 

150 

.05 

150 

.01 

300 

• 01 

300 

.05 


H+ 

CURRENT 

IN 

MICRO- A 

% 

DATA 122 6 

10.0 

-2.967 

10.0 

-6.940 

7.5 

2.507 

DATA 1230 

7.5 

0.759 

' 10.5 

-0.320 

1 1 .0 

-1.578 

1 1 .0 

-3.054 

11.0 

-4.815 

12.0 ‘ 

-5.210 

10.0 

-6.472 

10.0 

-4.430 

10.0 

-7.218 

DATA0 1 0 1 

8.0 

0.286 

8.0- 

4.766 

8.0 • 

2.462 

10.5 

0.520 

11.0 

-0. 1 53 

11.0 

-1.331 

10.5 

-2.836 

10.5 

-4.142 

1 1.0 

-5.237 

DATA0120 

11.5 

-3.619 

11.5 

-2.705 

10.0 

, 1 .493 

10.5 

0.614 

6.0 

-0 . 606 

6.0 

5.387 

5.5 

-0.084 

5.5 

-0.519 

5.5 

-2.988 

DATA0122 

6.5 

-0.838 

6.0 

0.800 

6.0 

0.937 

6.0 

I .952 

6.0 

0.857 

11.5 

1 .058 

11.5 

-1 .439 


S • E * 
OF 

MEAN 


0.693 

0.463 

0.649 


0.523 
0.605 
0.612 
0.439 
0.620 
0.679 
0.454 
1 .548 
0.240 


0.228 
0.918 
0.372 
0.578 
0.550 
0 • 564 
0.557 
0.587 
0.763 


0.278 
0.486 
1.086 
1.692 
1 . 959 
I .934 
0.148 
0.137 
2.534 


0.246 
0 . 1 30 
0.275 
0.846 
1 .568 
1.692 
0.91 1 



TABLE 13 COtITirUEO 


H+ 

HE 

H+ 


ENERGY 

PRESSURE 

CURRENT 

IN 

KEV 

MTGRR 

MICRO -A 

% 

300 

1.50 

11.5 

-1 .687 

300 

.01 

1 1 .5 

-2.373 



OATA0124 


450 

.01 

12.0 

-0.722 

450 

.05 

10.0 

-0.083 

450 

.60 

1 1 .0 

-0.274 

450 

1 .00 

1 1 .0 

-1.061 

4 50 

1.50 

1 1 .0 

-2.401 

450 

.01 

1 1 .0 

-2.336 

450 

.05 

1 0.0 

-1.158 

450 

1 .00 

10.0 

-3. 1 68 

450 

i .50 

10.0 

-4.224 

450 

.01 

12.0 

-2.661 

4 50 

• 05 

1 1 .0 

-3.547 


300 

.01 

300 

• 05 

300 

• 60 

300 

.01 

300 

‘.05 


300 

.20 

300 

.05 

300 

.01 

300 

.01 

300 

.05 

150 

• 05 

150 

.05 

150 

• 20 

150 

.01 

150 

.01 


450 

1 .00 

450 

.60 

450 

.05 

450 

.0 1 

300 

1.00 


DATA0220 

12.5 

3.270 

13.0 

0.158 

14.0 

-3. 103 

14.0 

2.700 

13.0 

-1 .673 

DATA04 1 4 

1 1 .5 

-2.180 

1 1 .0 

-2.381 

12.0 

2.213 

10.5 

-1.818 

1 1 .0 

-0.571 

6.5 

0.643 

6.5 

-0.085 

6.5 

1 .308 

6.5 

4.906 

6.5 

-0.239 

DATA0415 

11.5 

-6.377 

11.0 

-7.363 

11.0 

-2.289 

12.0 

-2.676 

DATA0513 

1 1 .0 

-1.1 12 


150 .01 


DAT*05!9 

^•5 - 6.200 


S .E. 

OF 

MEAN 

0.175 
1 . 3^0 


1.819 

0 . 7 70 

0.333 

0.310 

0.229 

1 .484 

1.049 

0.431 

0.1 52 

1.153 

0.603 


3.74 0 
1 .582 
0.227 
3.462 
1 .081 


0 • 669 

1.455 

2.585 

2.283 

1.345 

1.841 

2.194 

0.862 

3.769 

3.161 


0.172 
0.285 
2.189 
2 « 64 6 


0.228 


3.520 




TABLE 

13 CONTINUED 


H+ 

HE 

K+ 


ENERGY 

PRESSURE 

CURRENT 

IN 

KEV 

MTOER 

MICRO -A 

% 

150 

.01 

4.5 

0.667 

150 

• 01 

4.5 

-3.759 

300 

.20 

1 .0 

-3.706 

300 

.20 

1 .0 

-4.333 

300 

.20 

1 .0 

-5.943 



DATA0821 


300 

.20 

1 .0 

-0*994 

300 

.20 

1 .0 

-5.336 

300 

.20 

1 .0 

-2.316 

16 DATA PROGRAMS 

117 DATA RUNS 



S.E. 

OV 

MEAN 

2.725 
2 . 625 
2.275 
2.229 

1.7J2 


2.696 
2.346 
2.2 74 



TABLE 14 


RESULTS 

FOR EACH 

RUN, BY PAPA’-' 

JETER 


FROM 

H+ 



5 . E . 

DATA 

CURRENT 

IN 


OF 

PROG 

MICRO -A 

% 


MEAN 

! m 

K EV BEAM 

•20 MTORR 

HE 


DATA 1216 

4.5 

.2.760 


. 0.569 

DATA 1217 

5.0 

2.538 


0.900 

150 

KEV BEAM 

.01 MTORR 

HE 


DAT AO 1 20 

6.0 

-0.606 


1 . 959 

D..ATA0 J 20 

5.5 

-2.988 


2.534 

DATA01 22 

6.0 

0.857 


1 .568 

DATA0414 

6.5 

4.906 


3.769 

dataoau 

6.5 

-0.239 


3.161 

DA? AO 5 1 9 

4.5 

-6.200 


3.520 

DAT AO 5 1 9 

4.5 

0.667 


2 .725 

DATA05I9 

4.5 

-3.759 


2.625 

1 50 

KEV BEAM 

.05 MTORR 

HE 


DATA0 101 

8.0 

4.766 


0.91 8 

DATA0I20 

6.0 

5.387 


1 .934 

DA T A 0 1 2 2 

6.0 

1.952 


0 .846 

DAT AO A 1 4 

6.5 

0.643 


1 .84 1 

DATA 04 1 4 

6 .5 

-0.085 


2.194 

150 

KEV BEAM 

•20 MTORR 

HE 


DATA 12 16 

8.0 

1 .921 


0 .443 

DATA 1 2 1 7 

4.0 

2.624 


1 .205 

DATA 1225 

8.5 

1.791 


0.536 

DATA 1 226 

7.5 

2.507 


0.649 

data 1230 

7.5 

0.759 


0.523 

DATA0101 

8.0 

2.462 


0.372 

DATA04 1 4 

6. 5 

1.308 


0.862 



150 

KEV 

beam 

DAT A0 101 



8.0 

DATA0122 



6.0 


150 

KEV 

BEAM 

DATAO120 



5.5 

DATA0122 



6.0 


150 

KEV 

BEAM 

DATA0120 



5.5 

DATA0122 



6.5 


200 

KEV 

BEAM 

DATA1216 



8.0 

DATA 12 I 7 



7.5 

DATA 1 224 


10.0 


.60 

MTORR 

0.286 

0.937 

HE 

0.228 

0.275 

1 .00 

MTORR 

-0.084 

0.800 

HE 

0.148 

0.130 

1.50 

MTORR 

-0.519 

-0.838 

HE 

0.137 

0.246 

• 20 

MTORR 
1 .50 7 
0.675 
1.097 

HE 

0.495 

0,766; 

0.505 



TABLE 14 CONTINUED 


from h+ 

data current 

PROG MICRO-A 


DATA I 22 5 


i 1 .0 

DATA 1230 


10.5 

DATA01GI 


10.5 


250 

KEV BEAM 

DATA 1216 


8.5 

DATA I si. 7 


6.5 

DATA! 22 A 


18.0 

DATA 122 5 


10.5 

DATA 1230 


11.0 

DATA010I 


11.0 


300 

KEV BEAM 

DATA0I20 


10.5 

DATA0 1 22 


11.5 

DATA0 1 22 


11.5 

DAT A0 220 


12.5 

DATAG220 


14.0 

PATA0414 


12.0 

DATA04 14 


10.5 


300 

KEV BEAM 

DATA! 225 


10.5 

DAT A0 120 


10.0 

DAT A0 122 


11.5 

DATA0220 


13.0 

DATA0220 


13.0 

DATAP414 


1 1 .0 

DATA0414 


11 .0 


300 

KEV BEAM 

DAT A 1 2 1 6 


8.5 

DATA! 217 


8.0 

DATA 1224 


12V0 

DAT A 1225 


11.0 

DATA 1 22 5 


11.0 

DATA 1225 


19.0 

DATA1225 


4.5 

DATA! 225 


10.5 

DATA 1226 


10.0 

DATA 1230 


1 1.0 

DATA# I 0 1 


11.0 

DATA0414 


11.5 

DAT A0 51 9 


1.0 

DAT A0 5'1'9 


1.0 

DATA05I 9 


i .0 


S.E. 


IN 


OF 

% 


MEAN 

0.539 


0.374 

-0.320 


0.605 

0.520 


0.578 

•20 MTORR 

HE 


-0.308 


0.580 

-0.883 


0.795 

-0.785 


0.455 

-0.991 


0.430 

-1.578 


0.612 

-0.153 


0.550 

.01 MTORR 

HE 


0.614 


1 .692 

1 .058 


1 .692 

-2.373 


1 .340.. 

3.270 


3.740 

2.700 


3.462 

2.213 


2.585 

-1.818 


2.283 

.05 MTORR 

HE 


-1.321 


1 .240 

I .493 


1 .086 

-1.439 


0.911 

0.158 


1.582 

-1.673 


l .081 

-2.381 


1 .455 

-0.571 


U345 

•20 MTORR 

HE 


-1.359 


0.730 

-1 .279 


0.570 

-2.451 


0.527 

-3.909 


0.723 

-2.224 


0.451 

-3.221 


0.300 

-3.358 


0.699 

-3.850 


0.542 

-2.967 


0.693 

-3.054 


0.439 

-1 .331 


0.564 

-2.180 


0.669 

-3.706 


2.075 

-4.333 


2.009' 

-5.940 


1.712 
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■TABLE -14 CONTINUED • 


FROM H+ 

DATA current 

PROG MICRO -A 


DATA0821 


1 .0 

DAT A0 821 


1 .0 

DATA082I 


1 .0 


300 

KEV BEAM 

DATA 1 225 


11.0 

DATAQ120 


11.5 

DATA0220 


14.0 


300 

KEV BEAM 

DATA! 225 


1 1 .0 

DATA0I 20 


11.5 

DATA# 5 } 3 


1 1 .0 


300 

KEV BEAM 

DATA 122 5 


t 1 .0 

DATA0122 


11.5 


350 

KEV BEAM 

DATA 1 2 1 6 


9.5 

DATA 1 2 1 7 


9.5 

DATA1 224 


10.5 

DATA 122 5 


1 1 .0 

DATA 1230 


11.0 

DATA0101 


10.5 


400 

KEV BEAM 

DATA 1216 


9.5 

DATA1217 


9.5 

DATA] 224 


12.0 

DATA 122 5 


10.0 

DATA 1230 


12.0 

DATA01B1 


10. 5 


450 

KEV BEAM 

DATA# 124 


1 2.0 

DATA# 1 24 


11.0 

DAT A0 124 


12.0 

DATA0415 


12.0 


4 50 

KEV BEAM 

DATA 1230 


10.0 

DATA0124 


10.0 

DATA0124 


10.0 

DATA0124 


11 .0 

DATA041 5 


11.0 




S.E. 

IN 


OF 

% 


MEAN 

-0.994 


2.696 

-5*336 


2.346 

-2.316 


2.074 

60 MTORR 

HE 


-4.766 


0.175 

-2.705 


0.486 

-3.103 


0.227 

00 MT7RR 

HE 


-4.049 


0.302 

-3.619 


0.278 

-1.112 


0.228 

50 MTORR 

HE 


-2.896 


0.157 

-1.687 


0.175* 

20 MTORR 

HE 


-2.350 


0.407 

-3.142 


0.640 

-4.496 


0.669 

-3.304 


0.569 

-4.815 


0.620 

-2.836 


0.557 

20 MTORR 

HE 


-3.418 


0.672 

-2.194 


0.737 

-4.513 


0.529 

-4.585 


0 . 604 

-5.210 


0.679 

-A. 142 


0.587 

01 MTORR 

HE 


-0.722 


1.819 

-2.336 


1 .484 

-2.661 


1.153 

-2.676 


2.646 

05 MTORR 

HE 

f ' 

-4.430 


1.548 

-0.083 


0.770 

-1.158 


1.049 

-3.547 


0.603 

-2.289 


2.189 



TABLE 14 COim^UED 


from 

data 

PROG 


data l 2 1 6 
DATA 1 217 
DAT A 1 22 A 
DATA J 22 5 
DATA 1 22 6 
DATA! 230 
DATA01 01 


DATA 1230 
DATA0I24 
DATA0415 


DATA0 1 24 
DATA0124 
DATA04 1 5 


DATA0124 
DATA0 1 24 


H+ 


S.E. 

CURRENT 

IN 

OF 

MICRO-A 

% 

MEAN 


450 KEV BEAM 

.20 MT|RR 

HE 


9-5 

-3*815 


0.535 

9-5 

-3.477 


0.740 

l 1 -0 

-6.421 


0.696 

11.0 

-5.832 


0.886 

10-0 

-6.940 


0.463 

10*0 

-6.472 


0.454 

1 1 -0 

-5.237 


0.763 

450 KEV BEAM 

..60 HTORR 

HE 


10.0 

-7.218 


0.240 

11-0 

-0.274 


0.333 

11.0 

-7.363 


0.285 

450 'KEV BEAM 

1.00 MT^RR 

HE 


1 1 -0 

-1.061 


0.310 

10.0 

-3.168 


0.431 

11.5 

-6.377 


0.172 

450 KEV BEAM 

1 .50 MTORR 

HE 


11.0 

-2.401 


0.229 

10.0 

-4.224 


0.152 



TABLE 15 


7T Vs BEAM ENERGY 


AV. 


BEAM 

ENERGY 

(keV) 

BEAM 

CURRENT 

UA 

# OF 
RUNS 

# OF 

MEASUREMENTS 
OF it 

7T 

IN % 

S.E. 

OF MEAN OF tt 

tt FOR % CORRECTED 
FOR 

CASCADE AND TRAPPING 

100 

4.8 

2 

40 

2.65 

.53 

2.90 

150 

7.1 

7 

140 

1.91 

.27 

2.04 , 

200 

9.6 

6 

120 

.67 

.23 

.71 

250 

9.6 

6 

120 

- .78 

.24 

- .83 ■ 

300 

10.5 

10 

200 

-2.46 

.2 

-2.60 

350 . 

10.3 

6 

120 

-3.49 

.25 

-3.69 

400 

10.6 

6 

120 

-4.01 

.27 

-4.23 

450 

10.3 

7 

140 

-5.46 

.27 

-5.76 


Helium target gas pressure = 


2 mtorr 



Polarization Fraction 
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POLARIZATION Vs PROTON BEAM ENERGY 
(HELIUM TARGET GAS PRESSURE =.2mtorr) 



400 


Proton 

Beam 


Energy 

(keV) 



£ this work 
X Scharmann (1967) 
o Scharmann (1969) 


FIGURE 16 
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17, 18 and figure 17 show our results for tt (in %) as 
a function of He gas pressure at constant beam currents 
and at beam energies of 150, 300, and 450 keV respec- 
tively. The graph {fig. 17) uses the corrected it's.. 
Table 19 shows our results for r {in %) . as a function of 
beam current at 300 keV and .2 mtorr He pressure. 

5.2 Discussion of Results 

Let us first note from table 19 that a t-test of 
the significance of the difference of the mean tt*s for 
the normal and low beam current runs shows no significant 
difference. In other words the polarization is not a 
function of beam current. As can be seen from the table, 
this comparison was made at 300 keV beam energy and 
.2 mtorr He pressure. We chose to do the comparison at 
these values in order to obtain good statistics. The 
beam current was lowered by a factor of ten in order to 
make the comparison. The Van de Graaff accelerator is 
not capable of sustained runs at higher currents, nor is 
the Faraday cup capable of absorbing higher current runs, 
therefore we were not able to compare our normal current 
run with other than just the one 19.5yA run of 
December 25, 1972. A t-test between the normal and high 
current runs also implies no significant difference in 
means. However we must realize that there was only one 
19yA H + run and we should be leery of this single result. 



TABLE 16 


tt Vs HELIUM PRESSURE 

AV. 

He BEAM # OF tt FOR % CORRECTED 

PRESSURE CURRENT # OF MEASUREMENTS 7 r . S.E. FOR 1 

mtorr \iA RUNS OF tt IN % OF MEAN OF tt CASCADE AND TRAPPING 


.01 

5.5 

8 

160 

- .92 

1. 

- .97 

.05 

6.6 

5 

100 

2.53 

.75 

2.67 

.2 

7.1 

7 

140 

1.91 

.27 

2.04 

.6 

7. 

2 

40 

.61 

.8 

.66 

1.0 

5.8 

2 

40 

.36 

.12 

.39 

1.5 

6. 

2 

40 

- .68 

.14 

- .74 



Proton 

beam energy 

=* 150 keV 




& 


m 



He 

PRESSURE 

mtorr 

AV. 

BEAN 

CURRENT 

\iA 

# OF 
RUNS 

TABLE 17 

7T Vs HELIUM PRESSURE 
# OF 

MEASUREMENTS ir S.E. 

OF tt IN % OF MEAN OF it 

7r FOR % CORRECTED 
FOR 

CASCADE AND TRAPPING 

.01 

11.8 

7 

140 

. 8i 

.96 

.84 

.05 

11.4 

7 

140 

- .82 

.48 

• - .86 

.2 

10.5 

10 

200 

-2.46 

.2 

-2.6 

.6 

12.2 

3 

60 

-3.53 

.22 

-3.76 

1.0 

11.2 

3 

60 

-2.93 

.23 

-3.13 

1.5 

11.3 

2 

40 

-2.29 

.15 

-2.46 


Proton beam energy 


300 keV 



TABLE 18 


7 r Vs HELIUM PRESSURE 


AV. 


He 

PRESSURE 

xntorr 

BEAM 

CURRENT 

0A 

# OF 
RUNS 

# OF 

MEASUREMENTS 
OF TT 

tt 

IN % 

S.E. . 

TT IN % CORRECTED 
FOR 

CASCADE AMD TRAPPING 

.01 

11.8 

4 

80 

-2.1 

.92 

-2.19 

.05 

10.4 

5 

100 

-2.3 

.62 

-2.41 

.2 

10.3 

7 

140 

-5.46 

.27 

-5.76 

.6 

10.7 

3 

60 

-4.95 

.46 

-5.29 

1.0 

10.8 

3 

60 

-3.54 

.34 

-3.79 

1.5 

10.5 

2 

40 

-3.31 - 

.2 

-3.55 


Proton beam energy 


450 keV 



M1WL 


ill» 



1.2 

□ 

5 


300 keV 


— i — ♦ He 
*.16 Presst 
mtoi 


a 

5 


r 450 keV 



.2 

Q 


* He 
Pressi 
mtor 


.15 


TABLE 19 

TT Vs BEAM CURRENT 


AV. 

BEAM 

CURRENT 

yA 


# OF 

# OF MEASUREMENTS ir S.E. 

RUNS OF TT IN % 


19.0 

1 

20 

-3.22 

.3 

10.5 

10 

200 

-2.46 

.2 

1.0 

6 

120 

-3.77 

.88 


t-test Comparisons 


Low Vs . Normal 
High Vs . Normal 


t = 1.80 = 318 

t = 1.53 V = 218 


TT IN % CORRECTED 
FOR 

CASCADE AND TRAPPING 

-3.4 

- 2.6 

-3.98 


108 
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In any case we have shown that the polarization of the 
light from the 3 1 P-2 1 S transition in proton beam excited 
He is independent of beam current (over the range ^1 
to ^20. ii A H + ) . We are now able to compare measurements 
of tt . made at different beam energies and unavoidable dif- 
ferent beam currents. The accelerator is most efficient 
at 300 keV, and less efficient at 150 and 450' keV. We 
were therefore unable to run at the optimum of >10pA for 
all beam energies. The accelerator maximum currents are 
shown in table 20# Ten micro amperes was chosen as an 
optimum for several reasons? signal to noise ratios were 
better with this high current; . and the Faraday cup 
could tolerate this current for sustained periods of time. 
When possible we ran at M.0yA H+ , otherwise (at 100, and 
150 keV) we ran with as much current as we could get. 

We may now compare polarization measurements made 
at different beam currents having shown that the data 
taken at these different beam currents are compatible. 

The energy dependent results have been shown in table 15 
and in figure 16 (Corrected it Vs Beam Energy). The large 
error on the 100 keV point is due to the fact that we 
were only able to get two runs at that energy and pres- 
sure. The Van de Graaff accelerator, upon the installa- 
tion of a new (and slightly larger) exit canal at the 
source bottle, ceased to operate at 100 keV. it was 
deemed impractical to attempt modification of the ac- 
celerator in order to refine measurements of that one 
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TABLE 20 

APPROXIMATE BEAM CURRENT MAXIMUMS 


BEAM i BEAM 

. ENERGY ’ CURRENT 

keV p-A . 


100 

150 

200 

250 

300 

350 

400 

450 


5 


8.5 


12 

13 

20 

15 

13 

12 
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point. The points in figure 16 generally form a smooth 
gently curving line with the exception of the 400 keV 
point. Theoretical predictions (Appendix B, Table B2) 
from the work of Bell (19 61) are, except for the distor- 
tion. approximation at 178 keV, in very poor agreement 
with our results. Previous experimental work however is 
in good agreement with ours. We have abstracted the re- 
sults of Scharmann (1967), and Scharmann (1969), from 
their graphical presentations (neither of these works 
either tabularized their results or included error bars) 
and reported the estimates in table 21 along with the 
tabularized results of Van den Bos (1968). In the small 
region of overlap (100-150 keV) we find ourselves in 
much better agreement with Scharmann (especially con- 
sidering the curve shapes which we have not shown). 

Van den Bos* polarization measurements were done using a 
monochrome ter to select the X 5016 A line and he may not 
have completely corrected for the polarizing influence 
of the monochrometer. He also may have not taken care 
with regard to He pressure. He implies that he used 
.2 mtorr but does not explicitly say so. 

All other previous experimental work on the 
polarization of the X 5016 A line has been done at pres- 
sures .2 mtorr. This work is the first to extend that 
measurement to lower pressures, see summary in table 22, 
It has been difficult to do so since the signal to noise 
ratios fall off rapidly as one lowers the He target gas 





TABLE 

21 





COMPARISON WITH OTHER EXPERIMENTS 


i ■ 

BEAM 

ENERGY 

(keV) 

• THIS WORK 
. 2 mtorr 

VAN DEN BOS 
NO PRESSURE 
STATED 

SCHARMANN 

(1967) 

. 2 mtorr 

SCHARMANN 

(1969) 

. 2 mtorr 

THEORY-BELL 

DISTORTION 

100 

2.9 % 

-1. % 

2.5 % 

5.5 % 

14 % 

150 

2.04 % 

4. % 



! . 

200 

• 71 % 




(178 keV) 2.6% 

250 

- .83 % 


-2 % 



300 

-2.6 % 




(316 keV) 13% 

350 

-3.69 % 




' 4 

400 

-4.23 % 


• 

-5.5 % 


450 

-5.76 % 


-6 % 

-8. % 

-22 % 


112 



113 


TABLE 22 

CORRECTED tt Vs PRESSURE 


PRESSURE 

mtorr 

150 

keV BEAM 

300 

keV BEAM 

450 

keV BEAM 

.01 

- . 97±1 

• 

. 84 ± . 96 

-2,191.92 

.05 

2 . 67± 

.75 

- .861.48 

-2.411.62 

.2 

2 . 04 ± 

.27 

-2.6 ±.2 

-5. 761. 27 

.6 

. 66± 

.8 

“3.761.22 

-5.291.46 

1.0 

.391 

.12 

-3. 131.23 

-3.791. 34 

1.5 

- .74 + 

.14 

-2.461.15 

-3.551.2 
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pressure. We consider our .05 mtorr measurements as suc- 
cessful, but have less faith in our .01 mtorr results. 
There is simply very little . light reaching our detector 
from the target chamber when the pressure is .01 mtorr. 

We are therefore not surprised to find the .01 mtorr 
points in f rgur e' iV T7 • to -be o f questionable value. The 
.05 mtorr and higher pressure points do lie on relatively 
smooth sensible curves and in fact a similar curve shape 
was reported by Scharmann (1969) . He graphically dis- 
played results for tt as a function of pressure over the 
range .2 mtorr to ^8. mtorr °( for beam energies of 100, 

450, 600, and 835 keV) . Each of his curves shows a 
minimum (a most negative it), and as energy increases, the 
minimum occurs at a lower pressure. We have observed the 
same phenomena for our 300 and 450 keV curves. Our 
150 keV curve does not show a minimum, however from 
Scharmann's curves, we would expect the 150 keV minimum 
to occur at ^2 mtorr, a pressure which we do not reach. 
Scharmann’s results also showed the polarization (at 
100 keV beam energy) changing sign as the pressure 
changed. . (Also for the 4 1 P-2 1 S line at 100 keV.) Our 
results show a change in the sign of tt (for 150 keV beam) 
at >1.2 mtorr, and imply (at 300 keV) a change in sign at 
>.015 mtorr. Scharmann does take notice of this change 
in sign of tt, however he offers no explanation and in 
fact states that he has no explanation for this phe- 
nomena. We discuss this problem in Appendix D. 
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This work is limited in a number of ways. The 
lowest pressure we are able to accurately measure is 
.01 mtorr, we would have liked to have used lower pres- 
sures. However we could just barely get data at 
.01 mtorr as so little light was received by our detec- 
tor. We are also limited by the heat dissipation ca- 
pability of our Faraday cup, our Signal/Noise ratio 
could have been improved by using higher beam currents. 
(At some energies the beam current is limited by the 
Van de Graaff accelerator rather than the Faraday cup.) 

Improvements in this work could be made by . , , 

1. Obtain a new pressure measuring system 
capable of accurately measuring pressures 
to 10" 3 or 10~ 4 mtorr. 

2. Rebuild the Van de Graaff beam tube vacuum 
system so that one could have a differen- 
tially pumped target chamber at 10" 3 or 10" 4 
mtorr while the beam tube was held at 10" 5 
mtorr. (Now the beam tube ,is generally held 
at 10" 3 mtorr. ) 

3. Rebuild the target chamber electrostatic 
focusing system and entrance slits such that 
we get a beam of small cross section through 
the axis of the chamber, and such that the 
detector can be brought much closer to the 
beam region of the chamber. 

4. Modify the accelerator as so to be able to 
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obtain much higher beam currents especially 
at low energies. 

5. Devise a way to cool the Faraday cup while 
still keeping if electrically insulated. 

6. Use. the PDP-8 or equivalent to automaet the 
data taking. 

7. Shield the target chamber with mu metal so as 
to determine whether tt is in fact effected by 

• stray magnetic fields. Note tt should not be 
so effected. 

5.3 Conclusions 

It is clear from this work that the polarization 
of the light from the 3^P-2^S transition in proton beam 
excited Helium is still pressure dependent even at pres- 
sures as low as .05 or .01 mtorr. The results of 
Scharmann, as far as they go, bear this out. We are 
faced with a phenomena which exhibits poor- agreement be- 
tween theory and experiment, however we conclude that 
even we have not yet reached low enough pressures to find 
the natural value of the polarization. It remains for the 
next worker in the field to extend the measurement ' to yet 
lower pressures in search of a pressure independent value 
for the linear polarization fraction. 
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APPENDIX A 

DESCRIPTION AND MEASUREMENT OF POLARIZED LIGHT 
A.l Polarized Radiation 

It is well known that electromagnetic radiation 
behaves as a transverse wave. This means that the 
directions of both the electronic and magnetic field 
vectors of the wave are at all times perpendicular to the 
direction of propagation of the wave, as well as perpen- 
dicular to each other. Therefore, in order to completely 
specify the nature of a wave, one must among the many 
parameters , ' include those which describe the direction of 
either the electric or magnetic field vector as a function 
of time and position. Conventionally, one specifies the 
direction of electric field vector. We will from this 
point on use the term "light" in place of "electromagnetic 
radiation." Of course it is understood that light (in- 
cluding infrared, visable, and ultraviolet) is the name 
applied to a very small section of the electromagnetic 
spectrum. Furthermore comments made about polarized 
light do usually apply to other electromagnetic radia- 
tions . 

Light can be described as either unpolarized or 
polarized. Unpolarized light has no preferred direction 
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for its electric field vector, nor any preferred rota- 
tion of its electric field vector about the propagation 
direction. Conversely, polarized light does have pre- 
ferred direction and/or rotation for its electric field 
vector. As a matter of convenience, one can discuss 
three types of polarized light: 1. linearly polarized, 

2. circularly polarized, and 3. elliptically polarized/ 

1. Linearly polarized light. When the electric 
field vector of a light beam maintains a fixed direction 
in space, then the beam is called linearly polarized. 

See Fig. A.l. This is by far the simplest case of polari- 
zation and in fact is a type of polarization often produced 
in nature. 

2. Circularly polarized light. Light is said to 
be circularly polarized when the magnitude of the electric 
field vector appears constant while the vector itself 
rotates about. the propagation direction. There are of 
course 2 types of circularly polarized radiation. The 
conventional distinction is that if the locus of the "tip 1 * 
of the electric field vector describes a right-handed 
helix (such as the thread of a typical machine screw), 
then one has a wave of positive helicity (in optics 
called a left circularly polarized wave). Conversely, 
right circularly polarized light has the "tip” of its 
electric field vector describing a left-handed helix ; 

and possesses negative helicity. 

3. Elliptically polarized light. This is the 
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LINEARLY POLARIZED ELECTROMAGNETIC WAVE 



FIGURE Al 
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more complicated of the three cases due to the fact that 
both linear and circular components are present. To 
describe elliptically polarized radiation one must specify 
the handedness and magnitude of the circularly polarized 
component as well as the direction and magnitude of the 
linearly polarized component. Note that the major axis 
of the ellipse lies in the same direction as that of the 
linearly polarized component. An alternate equivalent 
description of elliptically polarized light is that it 
is made up to two linearly polarized components with 
differing amplitudes and non-0°, 90°, 180° phase angle 
differences. One must then specify each amplitude as 
well as the phase angle difference. 

For purposes of description one can draw sectional 
patterns for a beam of light. These drawings represent 
the electric field vector as seen by an observer located 
on the z axis (assuming that the beam is traveling in the 
+z direction) and looking toward the source of the radia- 
tion. See fig. A. 2. We might add at this point that a 
monochromatic beam has been assumed for the sake of 
simplicity. 

A. 2 Mathematical Descriptions of 
Polarized radiation " 

There are numerous equivalent mathematical descrip- 
tions of monochromatic waves, e.g. 
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Je| = ( S'/v(<a>* - ^ f A] 

or in complex notation 

!£/=;&[ 

or 

lef = P e .[e^e"*J 

where | E J is the magnitude of the electric field vector , 
w is the 2tt times the frequency of the wave, z is the 
position along the axis of propagation, A is the wave- 
length, and t is the time. The last equation is just an 
abbreviated form of the second where 

0 r 

Let us also state that Ee 1 ^ is called the complex ampli- 
tude, <p is called the phase angle (at t and z) and that 

the intensity of the wave is proportional to (E ) 2 . 

max' 

The direction of the electric vector has not yet been 
specified. Therefore the polarization of the wave has 
not yet been specified. 

At first glance one may easily arrange a system 
to describe a linearly polarized beam by simply choosing 
the x axis along the direction of the beam’s electric 
field vector. Then using the unit vector £ 

E = S'- (wf - f 

9 = Re [£»., a 
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l . 

One can then in turn describe circularly polarized 
light as the addition of two components of linearly polar- 
ized light of equal E max # directed along the x and y axis 
respectively, with a phase angle difference of +90°. By 
convention if the difference in phase angles, called y 
(y - #y“$ x ) / +90° the light is lef t-circularly polar- 

ized. Conversely, if y = -90° the light is right- 
circularly polarized. 


When the x and y linearly polarized components 
have different amplitudes, and when -180° < y < 180° 

(Y f 0), the beam exhibits elliptical polarization. If 
-180° < y < 0° then the handedness is right, conversely 
the wave is left-elliptically polarized if 0° < y < 180°. 

If one considers elliptical polarization to be 
the general case, then linear polarization results when 
Y = 0°, and circular polarization results when [y\ = 90° 


and (E ) = (E ) . 

max' x ' max y 

Strangely, it is very difficult to properly de- 
scribe a beam of unpolarized light. One can demand that 
no long term preference exist for handedness or linearity. 
(This precludes a "pure" coherent monochromatic wave.) 


Perhaps it is best to describe unpolarized light in terms 


of what it will and will not do (i.e., and operational 
definition) . When a beam of unpolarized light is divided 
into two completely polarized components (by an ideal 
analyzer) the components will deliver equal power. 

In the real world light is seldom either completely 
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polarized or unpolarized. One roust in such a case specify 
the degree of polarization. This is usually accomplished 
by dividing the beam of light into orthogonal polarization 
components (e.g. right and left circular; horizontal and 
vertical linear, etc.), measuring the intensity of each 
of these components, and then calculating the degree of 
polarization V from 



where 1^ and 1^ are the maximum and minimum intensities 
of the orthogonally polarized components respectively. 

Note that if one is discussing linearly polarized light, 

V reduces to tt the linear polarization fraction. 

.A number of different mathematical models are in 
use for the treatment of beams of polarized light. . These 
methods range from the geometrical methods of the Poincare 
sphere to the matrix calculus of the Stokes vector and 
Jones vector. 

The Poincare Sphere model uses a mapping technique 
to describe the polarization form of a completely polar- 
ized beam. It is most useful for predicting the change 
in polarization form as the beam passes through various 
active polarizing devices (analyzers, retarders, etc.). 

In this model , every point on the surface of a unit sphere 
represents a fferent and specific polarization form. 

By definition -he "north "pole represents lef t-circularly 
polarized light, and the "south" pole right. All points 



on the equator represent linearly polarized light with 
various inclination angles (to the equator) . One loca- 
tion on the equator (the prime meridian) is given the 
designation H for linearly polarized light parallel to 
the equator. Half way around the sphere (on the equator) 
the light is linearly polarized perpendicular to the 
equator, and given the designation V. At any location 
away from the poles and equator, one has elliptically 
polarized light. See Fig. A. 3. To specify the form of 
a completely polarized beam, one need only specify the 
angles X and ca (where 2 A = longitude and 2to = latitude) . 
This, is easily done since A is the azimuthal angle of 
the polarization ellipse and tan|cu| is the ellipticity 
(b/a). See Fig. A. 4. Furthermore the handedness is 
defined as left for the northern hemisphere (u><0) and 
right for the southern hemisphere (oo>0) . ■ 

When dealing with completely polarized beams and 
various retarders one can use this Poincare sphere to 
predict the change in the polarization form of the beam. 
This model is limited however in that one must be dealing 
with completely polarized beams. 

The Stokes Vector model of polarized light is 
perhaps the most versitle and useful. The Stokes vector 
is a four parameter column vector 



the poincare" sphere 


n 



FIGURE A3 
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ELLIPTICALLY POLARIZED LIGHT 


Y 



FIGURE A4 
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where all four entries have the units of intensity. 
Each of these parameters is defined in terms of a. „ 


measurement made with a set of ideal polarizing filters. 
This set of four ideal filters have the following proper- 
ties. All filters have a transmittance of .5 for un- 
polarized light. Then ... F-^, is a neutral isotropic 
filter; F 2 is a linearly (horizontally) polarizing filter? 
F 3 is a linearly (+45°) polarizing filter? and F ^ is a 
right-circular ly polarizing filter. A polarization • form 
independent detector is placed behind each of the filters 
in succession and four measurements are taken. The re- 
sult of each measurement is multiplied by l/.5.and the 
four values V^, V 2 , V 3 , and V 4 are reported. The stokes 
column vector is then calculated as follows 




An alternate and equivalent approach is to use 
the parameters from the wave description of light. One 
may then write (note time average ^ 

i = + ( f 4) 

m = 

s = <2 ijf> 


It is also conventional to normalize the Stokes vectors. 
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This is accomplished by dividing all four entries by I 
(which is the largest). The Stokes vectors are very use- 
ful in that they can be applied to partially polarized 
beams, monochromatic as well as polychromatic beams, and 
incoherent beams . Furthermore, the Stokes vector can be 
manipulated (by Mueller calculus to be discussed later) 
so as to predict the behavior of a beam when passing 
through retarders and polarizers. 

The Jones Vector model is applicable to situations 
in which the coherence and phase of the beam is important. 
For a specific location along the propagation direction 
(the z axis) the Jones vector is written as a two -element 
column vector 


r 

* for V- UJ t ) 


p) 

e 


( y 

* (fy . + t ) 

— 

n 

e 

_ 





where t and e are the phases of the x and y components 

x y 

at t - 0. There is a normalized form of the Jones vector 

■ 2 2 

where m + n =1. Conversion of this vector to a 
recognizable form is a bit cumbersome and in general use 
of the Jones vector is confined to problems in which 
phase is important. 

. Table A.l shows various polarization forms and 
the corresponding Stokes and Jones vectors. 

The Mueller Calc- us methodology provides a 
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means for manipulating the Stokes vectors so as to pre- 
dict the behavior of a_ beam of light as it passes through 
any series of retarders, and or polarizers. The method 
itself is of a matrix-algebraic nature. For this work 
it is sufficient to say that every retarder and or 
polarizer at whatever orientation can be represented by 
a 4 x 4 Mueller matrix. To determine the polarization 
form and intensity of a beam exiting from an " active” 
filter one multiplies the Stokes Vector for the incoming 
beam by the Mueller matrix of the "active" filter. The 
new Stokes vector then describes the exit beam. Using 
the methods of matrix algebra one can easily predict the 
result of any combination of "active" filters. 

The Jones Calculus methodology is similar to , the 
above except that a smaller (2 x 2) matrix is used. This 
smaller matrix is usually made up of complex elements and 
can preserve phase information. 

Standardized Mueller and Jones matrices are 
available in many works, notibly in Schurclif f (1966) . 

A. 3 Polarized Light and Matter 

When a beam of light (polarized or not) interacts 
with matter, the result of the interaction depends upon 
both the polarization of the initial beam, and the 
specific nature of the matter involved. 

Let us first consider (classically) the scattering 
of an unpolarized light wave incident upon a molecule. 
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The electric charges in the molecule will try to "follow” 
-the incident electric fields. As a result/ the molecule 
will radiate, i.e. scatter the light. The polarization 
of the scattered light will depend upon the location of 
the observer. If the incoming light is linearly polar- 
ized then the light is scattered anisotropically , and 
again the polarization depends upon the location of the 
observer. See fig. A. 5. 

Polarization can also be caused by reflection and 
refraction. Again the electric field of the incoming 
light will cause the charges within the material to os- 
cillate, thus multiple scattering is observed. Consider 
an unpolarized beam of light % (in air) incident upon a 
more optically dense material (glass). Applying the 
bo.undry conditions that the tangential electric and 
magnetic fields at the surface are continuous, one arrives 
at the situation shown in fig. A. 6. The degree of polari- 
zation depends upon the angle of incidence 0^. When 
+ 9 r - j the polarization of the reflected and re- 
fracted components is almost total. The incidence angle 

for which this occurs is named Brewster's angle 0 O , de- 

n 3 

fined by Tan 0^ = — . Where n 2 and are the indices 

of refraction for the two media. 

Certain materials exhibit the phenomena of 
Birefringence (double refraction) when interacting with 
light. Some crystalline solids fall into this category, 
that is they are optically anisotropic. In such a crystal 
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POLARIZATION BY REFLECTION AND REFRACTION 



FIGURE A6 
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(calcite or quartz for example) an incident ray of un~ 
polarized light will have two refracted rays (hence 
double refraction) as well as the usual reflected ray. 

One of the refracted rays will obey Snells Law. . . 

P, By = P+ Sr 

This ray is then called the ordinary or 0 ray, while the 
refracted component which does not obey Snells Law is 
called the extraordinary or E ray. See fig. A. 7. In 
1678 Huygens- discovered that the E and 0 rays were 
linearly polarized, and orthogonal. Birefringence is 
basically due to the fact that the velocity of the 0 ray 
is isotropic and therefore the material has the single 
index of refraction n Q for the 0 ray. The velocity of 
the E ray depends upon direction and therefore so does 
the index of refraction n for the E ray; The velocity 
of the rays depends upon the orientation of the electric 
field (the polarization) of the ray. This difference in 
velocity between the 0 and E rays enables the experimenter 
to selectively retard one of the linear polarization 
components. (This use of a birefringent crystal requires 
that the crystal be cut in a specific way... see Jenkins 
(1950) - A crystal arranged in this manner is called a 
retarder and has the properties of resolving a beam into 
two components, transmitting the components at different 
velocities, and recombining them with a phase difference 
that depends upon, among other parameters, the .thickness 
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of the retarder. One can then generate elliptically 
polarized. light by sending a beam of linearly polarized 
light through a retarder whose fast direction is at some 
angle 0 with respect to the initial plane of polarization, 
and whose retardance (phase angle difference) is 6. The 
results will be elliptically polarized unless G = + 45°, 

<$ = + 90° (circular); or 0 = + 90, 6 = + 90 (linear) or 
if 6 - 0°, 180° (linear). A retarder then is just a 
polarization form converter. 

The most important interactions of light with 
matter (for the practical use of polarized light) are 
dichroic interactions . Materials which exhibit dichroism 
will absorb one polarization form and transmit the 
orthogonal form, Polaroid sunglasses are a typical ex- 
ample of dichroic material.. The lenses are oriented so 
as to pass linear ly-vertically polarized light, therefore 
the horizontally polarized components are absorbed. But 
glare from the surface of a road, or body of water, or 
from many other horizontal reflecting surface is hori- 
zontally polarized, and removed by the dichroism of the 
sunglass lens. „ Although there are many materials which 
can be made to exhibit dichroism, the most commonly 
available is Polaroid Corporation's type H. This material 
is manufactured by alligning iodine molecules in a stretched 
sheet of polyvinyl alcohol. When a beam of light passes 
through the sheet, the component of the beam whose electric 
field lies along the long axis of the iodine molecules 
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will be absorbed. The orthogonal component is trans- 
mitted (with some attenuation due simply to the optical 
density of the sheet) . Land (1951) has done a review of 
the history of dichroic polarizers (He is the inventor 
of the modern dichroic polarizer) . 

Optically active materials such as Turpentine# 
Quartz crystals and sugars (to be optically active a 
sample of the compound must contain a significant majority 
of one of the mirror image versions of itself; either the 
right or left handed version) can rotate the plane of 
polarization of a beam of linearly polarized light. 

Jenkins (1950) describes- Fresnel ' s explanation for 
optical rotations. Fresnel assumed that optically active 
material treats the entering linearly polarized light as 
tv/o equal amplitude# zero phase difference, counter- 
rotating circularly polarized components. He further 
assumed that the propagation speed was different for the 
Right and Left components; producing a non-zero phase 
difference. Then because of the non-zero phase difference# 
when the beam exits the material and the two counter- 
rotating components recombine#, we once again have linearly 
polarized light# but the plane of polarization has been 
rotated by an angle 6/2 (6 is the phase difference) . It 
is apparent then that the magnitude of the rotation 
depends upon the amount of optically-active material the 
beam passes through. Although Fresnel's explanation does 
agree with experimental evidence it does not explain 
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the process on the molecular level. That explanation 
was first put . forth by Reusch; he suggested that the 
molecules in an. optically active substance are aligned 
in a helical pattern. This was later verified in other 
experiments . 

A. 4 Semiclassical and Quantum Mechanical 
Aspects of Polarized" Light ~~~ 

In the previous section we have discussed the 
interaction of light with matter, without regard to how 
the light was created. We will now rectify that omission 
by discussing (both semi-classically and quantum mechani- 
cally) the origins of light within the most elementary 
radiating systems. 

We begin with an electric dipole. Consider a 
fixed positive charge with a negative charge in linear 
periodic motion about it. Let us assume that our observa- 
tion point is located a distance r Q from the positive 
charge with r^ ~ r Q >> r. See fig. A. 8. We can then . 
write the retarded scalar and vector potentials (for the 
observer) of the moving negative charge as follows. 



If we expand the above terms using 


A. 4.1 
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(l-x)~° 


I + nx + A - 1 +■ oix 7 ).-- 


2 ' 


> f 2 

and if we neglect terms of order (—) and higher... 

c 

then 


<P~ ~ 

' Tf > r « C \ t- f J. 


A. 4. 2 


^ r „ c (<-.4 


If we now neglect the first term in the scalar potential 

-Q. 

, o 

— -f-pi-rm f rnm f ha naci + 

r 


expression (at large r Q the ^ is cancelled out by the 
+Q ° 

— term from the fixed positive charge), and use the 


facts that 


B - - V d> - 1 ys 
r it 

it* 


' A. 4. 3 


we arrive at 


p - r 
* ^ * 

W -■ 7 ic *• £] 


A. 4. 4 


where is the component of the negative charges ac- 
celeration in a plane perpendicular to r Q . The observer 
(at sufficiently large r Q , in the wave zone) sees a plane 
polarized electromagnetic wave. Let us surround this 
oscillating and radiating di :-ie with an imaginary sphere 
(radius r Q ) such that the polar axis is along the dipole 
moment, the observed polarization direction of the wave 
will then depend upon where on the sphere the observer 
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is located. See fig. A. 9. Since the dipole makes an 
angle 6 with respect to r Q ., and since Qr = p * p cos(t)t)z 


. . . then 


and 




C 

Kc' 


A. 4. 5 


/e / - IHl - S * e A - 4 - 

where H is always parallel to the equator of our sphere 
and E is always perpendicular to said equator. Thus the 
wave is completely linearly polarized. 

Let us now consider the electric rotor where 
the negative charge travels about the fixed positive 
charge in a circular path. Rather than repeat the calcu- 
lations from the beginning, we can instead resolve the 
circular motion into two linear harmonic components with 
the same frequency and amplitude but with a phase dif- 
ference of +tt/2 • If we again use our imaginary sphere , 
this time with the polar axis along the symmetry axis 
of the rotor, again the form of the polarization will 
depend upon the location of the observer. For 0 = 0°, 
180° one observes a circularly polarized wave, at 9 = 90° 
one observes a linearly polarized wave, and anywhere in 
between one observes an elliptically polarized wave. 

There are of course higher order terms which we 
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have neglected, a summary of the radiation from higher 
electric multipoles as well as from magnetic radiators 
is available in Feof ilov (1961) . 

We can now discuss polarization in the normal 
Zeeman effect. Let us place our linear dipole oscillator 
in a magnetic field H not directed along the dipole 
moment. The equation of motion of the negative charge is 


ID r = - m u) B l r - | [ F * H J 


A. 4. 7 


which can be written as the three scalar equations (choos- 

ing the z direction to be along H) 

p) x r - m IAJ , 1 X - ? y M 

• n f 5 ~ i*» A. 4 .8 

mi* - - 

The z equation solves immediately the z = z cos w.t. We 

o o 

will try solutions of the form x = A e llot and y = A 

X J ' y t 

for the x and y equations. This procedure yields 
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therefore 


\aj ~ f H 

\J y/»* , C v “ 


A. 4. 11 
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However for magnetic fields , H , of realizable magnitude, 
-- w o ~ > ilc and theref ore. dropping, terms- of order (21L.) 
we arrive at 

u> = W, t AW 



A w> z 


<3 H 
z me 


A. 4. 12 


and note Aw << w . We may now solve for the A and A 

O • x y 

by again using x = A e la)t and y = K e ia)t in A. 4. 8 and 

a 2 y 

then removing terms of order .. We then find 


- r A Ay 


and the full x and y solutions are 

a {W t 1 4 Hi) t 


X = A* <2 


Y - e 


i [K, +A*>)t 7 


A. 4. 13 


A. 4. 14 


Therefore when w = w q + Aw, x leads y by a phase angle of 
tt/2 and the radiation emitted by this oscillator (when 
viewed from the z axis) is circularly polarized? and when 
w - u) Q -Aw, x lags y by a phase angle of tt/2 and the ; 
radiation is again circular butopposite in direction 
(called a components) . 

In summary, when a linear dipole oscillator is 
placed in a magnetic field one observes three frequencies 
of oscillation w q , w Q + Aw, and w q - Aw. Both the in- 
tensities and polarizations of these radiations depend ■ 
upon the location of the observer. If the observer is 
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on the z axis he will observe only circularly polarized 

components of frequency w + A w. (the linearly polarized 

radiation has zero intensity on the z axis). However, 

if the observer is on a line perpendicular to the z axis, 

in the plane containing the dipole, he will then observe 

all three radiations cu , w + Aw. Now however the a com- 

o o — 

ponents (w = w Q + Aw) will appear as linearly polarized 
with the electric field vector perpendicular to the mag- 
netic field direction. (Analogous to observing an electric 
rotor edge on) the w Q radiation (the it component) is of 
course linearly polarized with its electric field vector 
parallel to the magnetic field direction. 

This classical description of the normal Zeeman 
effect is limited to cases in which the source of radiation 
can be considered a simple system of linear dipole os- 
cillators. 

We now turn to quantum mechanics for a description 
of polarized light on an atomic scale. We will confine 
our discussion to electric dipole transitions where 
A A = + 1, and Am = 0, +1. These selection rules are 
necessary for non zero matrix elements between the upper 
and lower states. The transition probability is of course 
proportional to the square "of the matrix element. The 
selection rules can be demonstrated as follows. The matrix 
element 

d ** = 5 p ^ A - 4 - 15 
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(where p is the dipole moment of the atom electron system) 
oan.be resolved into the cartesian elements 

D* -- -e t *1- 

hi --- e- Sfl y t ^ a. 4.1 

d;„- 

assuming that the electron moves in a central force field, 
we can then write 




A. 4. 17 


where Jt is the orbital angular momentum quantum number, 
and m is the "z" component of £. There are 22, + 1 values 
of m, and m ranges by integers from -2. to +2,. Using 
A. 4. 17 in A. 4. 16 and with the z axis as the polar axis we 
arrive at . . . 
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The first two integrals will be zero unless 
(m f r m) » + 1 or Am « j* 1 . The last integral will be 
zero unless (m* — m) » 0 . Also note that when Am = + 1 
we have in effect linear dipole oscillators in the x and 
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y directions and circularly polarized light, (as viewed 
from a point on the z axis) these are again the o com- 
ponents. When Am - 0 we again. have the equivalent of a 
linear dipole oscillator with the E vector parallel to 
the z direction, which is the it component. Further note 
that when viewed from position in the xy plane one sees 
the a components as linearly polarized with E perpendicular 
to z and the tt component as linearly polarized with E 
parallel to z. It can also be shown that for Am = 0, + 1, 
A& must be + 1. These selection rules can be generalized 
to the set of quantum numbers rUjM. Where then Aj = + 1 , 

6 but j = 0 j = 0 is forbidden and AM = 0 , + 1. 

The above selection rules are appropriate in light 
of the following arguments... the angular momentum of 
linearly polarized light is zero, therefore when such 
light is. emitted, the change in the angular momentum com- 
ponent along any axis (mft) should not change, hence 
Am - 0. Circularly polarized light carries angular 
momentum + h one would therefore expect Am = + 1. 

We are now prepared to discuss the Zeeman effect 
quantum mechanically. Let us choose the z axis along 
the external magnetic field H. This problem has been 
solved in many books (see Bethe (1956) pg 208) with the 
following result... 

EntjH ~ A. 4. 19 

where - 
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and m takes on (2j 4- 1) values. The energy E' therefore 
takes on (2j + 1) values, one of which is still the zero 
magnetic field value ^ . The wavelengths, of the radia- 
tion are then 

y, a . t — - 

(M, f, H A.4.2 

The polarization rules are again, tt component when AM = 0 
(and therefore X - X Q ) , and a component when AM = + 1 
(and therefore X = X Q + AX) . The apparent polarization 
of the it and a components will depend upon the location 
of the observer. If the observer is located upon the z 
axis then the intensity of the it component will be zero 
and only the circularly polarized a components will appear 
(at wavelengths X - X Q + AX) . If the observer is located 

. ■ f 

in the xy plane, then both the 7r component (linearly 
polarized with the E vector parallel to z, with wavelength 
Xq) / and the a component (appearing a linearly polarized 
with the E vector perpendicular to the z axis at wave- 
lengths X = X Q + AX) will appear. As a specific example 
consider Helium with s = 0, therefore J = L and g = 1. 

See fig. A.10, It was also shown by Feojilov (1961) that 
for this specific , P 1 ^ ] s q transition the ratio of 
intensities for the it and a components (observed in the 
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NORMAL ZEEMAN EFFECT 



The presence or absence of it or a radiation 
is indicated schematically in the lower two 
lines. 


FIGURE A10 
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x, y plane) is 1:2:1 for o._^:ir:a + ^. 

The question now arises of what happens to the 
tt and a components as the magnetic field H goes to zero. 
Certainly the wavelength splitting decreases toward zero, 
and in fact the magnetic sublevels become degenerate. 
However (experimentally), the tt and a components do con- 
tinue to exist even at zero field. Heisenberg ' (1925) 
has stated the principle of spectroscopic stability... 
(paraphased) the state of polarization does not change 
when the external additional magnetic field, which is 
superimposed on the system in such a way that its symmetry 
remains unaltered, tends to zero. 

The entire above discussion has concerned itself 
with one elementary radiator (or one atom). In general 
for a macroscopic system of many radiators the polariza- 
tion will average to zero unless some anisotropy is intro- 
duced. In this work the anisotropy is introduced by the 
energetic proton beam. The beam direction is chosen as 
the z axis and observations are made from within the xy 
plane. See figure 2. 

A. 5 Measurement Methods for Polarized Light 

The devices which are used to analyze polarized 
light in fact apply those phenomena which we have dis- 
cussed in A. 3. Some of these devices are? the linear 
polarizer (dichroism) ? the polarizing prisms , Nicol , 
Glan-Thompson etc. (birefringence) ? retarders (bire- 
fringence) ? and the pile of plates polarizer (reflection 
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and refraction) . 

The linear polarizers are the. invention, of ,E. K. 
Land and generally the product of Polaroid Corporation 
(Cambridge, Massachusetts). The so-called H-sheet is the 
most common, and is covered by a number of patents. Dur- 
ing manufacture a thin sheet of polyvinal alcohol is 
heated, stretched and then for support purposes laminated 
to a sheet of cellulose acetate butyrate. The polyvinal 
alcohol face is then wet with an iodine rich solution which 
will leave a residue of aligned iodine molecules within 
the polyvinal alcohol. This sheet is then laminated be- 
tween plates of glass to form the completed linear polar- 
izer. These polarizers can be made in various colors and 
transmittances . In this work a type HN-32 implies a peak 
transmittance of 32% for unpolarized incident light. The 
behavior of this polarizer for normal incidence is cate- 
gorized in table A. 2. Note that ^ (the major principle 
transmittance) is defined as the ratio of transmitted to 
incident intensity when the linear polarizer is placed in 
a normally incident, linearly polarized beam of light 
oriented to maximize the transmittance. The minor princi- 
ple transmittance, K 2 , is obtained from the same ratio 
during minimum transmitance. The principle transmittance 

ratio R t is defined as k^k.,. For one polarizer R fc is on 
the order of 1.5 x 10* (at 5000 A). Reflection losses 
do occur at the surfaces of linear polarizers, however 
these losses (for normal incidence) are isotropic and 



TABLE A. 2 

BEHAVIOR OF THE HN-32 POLAROID ANALYZER 


Wavelength 

K r 

K 2 

3750 ■- 

.33 

.001 

4000 

.47 

.003 

4500 

.68 

.0005 

5000 

.75 

.00005 

5500 

.70 

.00002 

6000 

VO 

• 

.00002 

6500 

.70 

.00002 

7000 

.77 

.00003 

7500 

.84 

.0002 
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amount to typically <4%. 

There are of course many other dichroic linear 

polarizers which do not immediately concern us. De- 
tailed-information can be obtained from Polaroid Corpora- 
tion, Cambridge, Massachusetts. 

The linear polarizing prisms have been in large 
measure replaced by the dichroic sheet polarizers. How- 
ever one great advantage the birefringent polarizer has 
over the dichroic type is that a polarizing prism can 
separate and deliver both of the orthogonal linear com- 
ponents, there is no significant absorbtion. Very few of 
these prism polarizers however can be adjusted to leave 
the path of the incident beam undeflected, this is the main 
cause of their replacement by the linear dichroic polar- 
izers. 

Polarization form conversion is done with the use 
of retarders. The most common retarders are quarter and 
half-wave plates. The quarter-wave plate will convert a 
beam of linearly polarized light into a beam of circularly 
polarized light (when properly oriented) . The half-wave 
plate can be used to rotate the plane of polarization of 
linearly polarized light. There are also circular and 
elliptical retarders, so that in general with the proper 
combinations and orientations of linear polarizers and 
retarders one may convert from any one polarization form 
to another. 

The pile of plates polarizer is usually used for 
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production and/or analysis of vacuum ultraviolet polarized 
light. The^ polarizers previously mentioned generally are 
unusable in the ultraviolet region of the spectrum. There- 
fore the pile of plates reflection/refraction polarizer 
arrangement, even though extremely cumbersome, is grudging- 
ly used . 

Given a beam of light of unknown polarization form, 
one can use the devices previously discussed in order 
to analyze the polarization form. Once the form has been 
analyzed, specific polarizing analyzers can be used to 
separate the orthogonal components for the purpose of 
intensity measurements.' The results of these intensity 
measurements can be used to calculate the degree of 
polarization. The beam is then completely analyzed. In 
order to make the intensity measurements one needs a 
detector which meets two general requirements. Many 
different detector types can be used, the requirements 
being that the detector be sufficiently sensitive and 
also that the detector be completely polarization in- 
sensitive. 

In this research the HN-32 analyzer is used to 
isolate the linear polarization components so that we 
may use a photomultiplier detector system to measure 
the/ intensities of the components. We are then able to 
calculate the linear polarization fraction. 
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APPENDIX B 
THEORY 

B.l Collison Theory 

We must eventually use inelastic scattering 
theory to describe our process... 


• H + + -*■ H + + H* 

e e 

However, let us first review elastic scattering. 

We will assume the ideal case of a fixed scatter- 
ing center at the origin impacted upon by a beam of 
particles incident along the +z direction. We can re- 
present the scattering center by a potential V(f) , and 
let us assume that the incoming beam is monoenergetic . 

If after the elastic scattering occurs, the particles 
are detected far from the origin by a device which sub- 
tends an angle , then the differential scattering cross 
section is 


<nv) 


a/ j-a 


B.l.l 


where N is the flux of particles in the beam, and dN.i's 
the flux counted by the detector. The total scattering 
cross section g* is then 


c; - 


~ ^ J sl 


B.l. 2 
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We can describe the monoenergetic beam of parti- 
cles by a plane wave, with wave function - 

a <• z 


f = e 


B.1.3 


The wave function after the interaction will then be 

a k'i 


f = e 


-h 


V 


B.1.4 


where at large r, we must have 


V f (9, P) £ 


< K * t 


B.1.5 


i.e. an outgoing spherical wave. The particle flux can be 
calculated from the Wavefunction by using the expression 
for the probability current density... 


S = {?**) ?] 


B.1.6 


When this is applied to B.1.5 and we further assume large 
r , then 

tfw )! 1 


r . H ItML 1 c v - 


c v 


B.1.7 


ft V p 

where - - •*. V , the classical velocity of the parti- 
cle, and S is the scattered particle flux per unit area. 

Since dsi = <p and since an element of 

area in spherical coordinates is <1 A z c x s^ & J & <p we can 
then write 

cl A 


Iji * 


B.1.8 
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and therefore 




N F J A 


r r 


but N O' A -K- ' cl Af therefore 


Aj/ __ N r (e, <P) 
r 1 


B.1.9 


B.1.10 


This is the scattered particle flux per unit area at 
angles 6 , and <t > . 

We can apply B.1.6 to B.1.3 to arrive at the 
incident flux N 


Hi * 1* «ir 

' m 


using B.1.11 in B*1.1Q we get 


B.1.11 


la - or ir(6. p) B.1.12 

JA . r l 

a comparison of B.1.12 and B.1.7 yields the important 
result that the differential scattering cross section is 
equal to the square of the scattering amplitude, i.e. 


b. l .13 

The solution of many scattering problems is then 
reduced to finding the scattering amplitude f(&* fy . 

Recall the scattering amplitude is the coefficient 
of the spherical outgoing wave in the asymptotic form 
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+■ r A K*r 0 A K ' r 

f (r) - [e + *(*.?) t~ j 

when this is used in the Schrodinger equation 


B.l 


(yi- + k l ) vb J = ^ /(r-j tk) 


B.l 


one arrives at (see e.g. Merzbacher (1961)) 


VW) -- -J? VM ^(r)Jr' 


B.l 


As it stands, this result is not explicit and not useful 

. •+ 

The first Born approximation simply replaces y ( r ') by 
a plane wave £ k * r so that 


. 7 -' — ' 
-/i k >r 


ited = \ e V(f) e 


A k>r ' 


Jr' B.i 


This approximation is generally valid for short range 
weak potentials, or very high energy incident particle. 
The scattering amplitude is often written as 


¥*,*)* -^Je +vp V^) d r 


B.l 


where. a T , If/* Ik/ , ^ 


and 


os 0 x 


K*Z ' 

/«/* 


(see figure B.l) and therefore we have... 


[e %y v(r) Jr'l 


r 


B.l 


14 


15 


16 


.17 


.18 


.19 


and 
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(f t r ^ <T(Q,<P) 


-B.1.20 


Let us now go on to a discussion of inelastic 
scattering * The vectors X and K' no longer have equal 
magnitudes since energy is transferred to the target 
atom. The total system Hamiltonian is H where 


H = H + H + H 

atom particle interaction 


B.1.21 


In general, the transition amplitude from a state a to 
a state $ is given by 


<V> i H t i > B - 1 - 22 

For our situation, the ty's are product wave functions 




t -- 






B.1.23 


where o and n represent the ground and excited states. 
Conservation of energy requires 


rl («'-**) - U-e.) -- K -«>) 


B.1.24 


where E refers to the atom and W to the incident particle. 
The momentum transfer relations are now £ * •- ^ 

+ and , see 

again B.l. Equation B.1.22 now becomes, {in the Born 
approximation}... 



where A is representative of the atomic electron coordin- 
ates . 

To obtain the differential and finally the total 
cross sections for .the excitation, one must use B. 1.2-5 
in Fermi's golden rule... 




B.1.26 


The density of states factor, y> (^ n ) is equal to 
so 


(ini) 1 c v 



B.1.27 


and 




B.1.28 


At lower energies (or with target atom potentials 
of greater range) , the Born approximation becomes less 
accurate. One then turns to the Distorted Wave Approxi- 
mation. The initial and final wave functions of the 
particle are no longer the simple free particle wave 
functions . 
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In fact, the Distorted Wave Approximation will 
.use the _ approach of Eqs, J5 .1.2 2 and B.1.27, but with wave 
functions which have been distorted by the interaction 
potential- Details are available from Bethe (1968), 
from Mott (1965), and also from Hasted (1964). 

B.2 Cross Sections and Polarization 

The polarization of light emitted (in a particular 
direction) by an impact excited gas will depend upon the 
relative cross sections for excitation of the magnetic 
substates of the upper level. 

The methods of the preceding section can be 
generalized to calculate the cross sections for excitation 
to the various M levels. A conventional label for these 
cross sections is where a is the total orbital angular 
momentum quantum number L, and M represents the magnetic 

. Uft) . i 

substate. Let us next define f\ ^ (xj } ( rj , and 

("?) as the transition probabilities for emission of 
light from the Mth substate with the light’s electric 
field vector aligned with the x,y, or z directions re- 
spectively. Also, since our system is symmetric about 
the z axis 

i'Cm '/Cm * tl** - '/>*(*)] 

The intensity of light emitted with its electric field 
vector polarized parallel (perpendicular) to the z axis 



166 


is then 




, £ fl„ lY) Q„ 


X 


B.2.2 


7*1 


We can then write the linear polarization fraction 


.|r - A« fl (y)l q 

^ F/C + aZ (y)'] Q? 


0< 


B.2.3 


u. 




Using the Zeeman intensity formula from Condon (1967) for 
a J -*• J-l transition we get 

3,' *< (T l -M x ) 




B. 2 . 4 


where o 4 7T r anc/ (/) ^ ^*c 


then Eq. B.2.3 becomes (with use of B.2.4 and 7 ) 

for a 'p -* *5 transition 

7T r fi.2.5 

The theoretical problem has been reduced to finding 
the G?, and Q, for the *3 P state of Helium. 

Both the Born and Distorted Wave Approximations 
can be used to find the <£ and given a proper 

choice of wave functions, i.e., wavef unctions in which 
the M dependence has not been summed over. The matrix 



elements must be taken between various M states so that 
the transition probability will show any M dependence 
they might have. Mott (1965) gives an example of this 
type of procedure. Bell (1961) and others who have dohe 
this type calculation have used multi-parameter wave- 
functions whose parameters have in general been chosen 
to satisfy calculations of oscillator strengths . The 
cross sections arrived at by Bell (1961) are shown in 
table B1 as abstracted by this author from Bell's (1961) 
graphical results. Table B2 shows the values of tt as 
calculated from the cross sections of table B1 using 
equation B. 2 .5 



TABLE B1 

THEORETICAL EXCITATION CROSS SECTIONS 


Energy 
in keV 

Q o in 
Born 

7ra 2 
o 

D.W.A. 

Q ±i 

Born 

in Tra 2 
o 

D W. A 

100 

.0238 

.0198 

.0271 

■ / 

.0150 

178 

.0166 

.0175 

.0238 

.0166 

316 

.0107 

.0122 

.0192 

.0158 

450 

. 0080 

.0093 

.0163 

.0144 


! Cross sections for the excitation of He to the 

3 P 0 and 3 P ±1 states calculated by both the Born Ap- 
proximation and Distorted Wave Approximation (D.W.A.) 
methods. These numbers v;ere abstracted from the 
graphical results of Bell (1961). 



TABLE B2 


THEORETICAL LINEAR POLARIZATION FRACTIONS 


H + Energy 
in keV 

it Theoretical 
from Born 

tt Theoretical 
from Distortion 

100 

- .065 

+ .138 

178 

- .178 

+ .026 

316 

- .284 

- .129 

450 

- .342 

- .215 


from Bell's fraction * calculated 

rrom Bell s (1961) theoretical cross sections. 
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APPENDIX ' C 


co; iputfir. v no 


.1. la.; 


l-.io prog rain is used to calculate the behavio 
the interference filter. * 

LAM " ' 5,19 FRIDAY AUU 24, 1973 

.:^^.; , r N ;T.^ c f 8 F,LTER 

*20 PH] NT USING 130, N 
J30 t ANGLE LAMBDA 

U0 FOR A = 0 TO !0 STEP 2 
150 LET X-HALKA) 

nALrsir^;;:^"”" 8 ''"-- 8 ”- 


r of 


100 !#*0 
190 next a 
S00 PRINT 
810 NEXT ri 
280 END 





C . 2 \ T b:\TO OL 


This program is used to calculate it on a run by 
run basis according to equation 4.1.3. 


NtV/POL 9*22 FRIDAY AUG 24, JQ73 


1 00 DIM P(25) 

IIP! PRINT ’LINEAR POLARIZATION FRACTION’ 

120 PRINT •# OF RUNS TO ANALYZE IS': 

130 INPUT N0 

140 PRINT ‘LIST ALL P*i 

150 INPUT Ai 

160 FOR N 1=1 TO N0 

170 READ Ai,A2,A3,B0,Bl 

180 LET C / D=-0 

190 FOR N2 = 1 TO B0 

200 -READ B2,B3 

210 LET C=C+B2/B0 

220 LET D=D+B3/60 

230 NEXT N2 

240 IF AS = ' YES * THEN 270 
250 LET Z = 2 
260 00 TO 310 

270 Let z = i 

280 PRINT 
290 PRINT 

300 PRINT • V C H C 

310 LET F = 0. 

320 FOR N3 = l TO B1 
330 READ V , N , H , M 
340 LET A=(H-C )/(M-D> 

350 LET B= (V-C )/ <N-D> 

360 LET P<N3)=<A*B)/(A+B> 

370 LET K=F+P(N3)/B1 
380 0 0 TO 3 90,4 10 ON 7. 

390 PRINT USING 400 , V, N, H, M, P < N3 ) 

400 t 9*0f* 0*000 00900 00000 09.0000* 

410 NEXT N3 

420 LET S = 0 

430 FOR N4 = 1 TO Bl 

440 LET $=S+ CP<N4 )-F> t2 

450 NEXT N4 

460 IF Bl =1 THEN 490 

470 LET S=SQR(S/CB! -1 ) > 

480 00 TO 500 

490 LET S =0 

500 LET E = S/SURCB1 ) 

510. LET 1=0 

520 FOR N 5= 1 TO Bl 

530 LET 0 =ABS <P <N5 ) -F ) 

540 IF G < 2 * S THEN 560 


P* 
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550 LET I •- } + j 

■560 NEXT l,-r> 

570 PRINT 

580 PRINT I/S IN.(j 5 90 > A ! * A? > A 3 j B 1 

590 i*p* KEV, 0mff1f □ H£> A H+ 

60P PRINT USING 610, S, I, E 

610 t ST DV -***#### DM > P <;ir 

620. LET F0=F-E. 2 SIC 

630 LET F 1 = F + E 

640 PR I NT USING 650, F0,F,F1 

650 J <P-SE> s:0t# NEAN P =#* 

670 NEXT N* 

680 EN» 


*» trials 

S *E* =.# ,0 tint* 

< P+ SE )=#,#*<## <? 
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■ C» 3 ST AT 

Thi-j program is used to calculate tt f or ' a 
collection pi runs ( from results stored in a ci le 


STAT 


1 00 
110 
1 20 
130 
140 
i 50 
1 60 
170 
180 
, 193 
200 
210 
220 
230 
24 0 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
4J fl 
420 
430 


9:25 FRIDAY AUu 24/ 1973 
DIM PC300) 

u^put At 11- * T ° BE AliALY7ED Ifi,; 

OPEN 1 /AS/ INPUT 
GET f-!/(g 
LET F s=0 
For Jsrj TO M 
GET P< I > 

LET F*F+p<i)/p 

NEXT I 

LET S - 0 

FOR 0=1 TO M 

LET S =S+ (P( J)*F) r 2 

NEXT J . 

LET S=S0R(S/(M-J )) 

LET K~0 

for L ~ 1 TO M 

LET G=ABS ( P (L ) -F ) 

IK u<2*s then 300 

LET K =K + 1 
NEXT L 

LET E=S/SOR<m> 

LET F0 =F -E 

LET F I s F + E 

PRINT USING 350/M 

\*** P’S IN CALCULATION 

PRINT USING 370/S/K/E 

J ST DV =*.###** ### pTs 

PRINT using 3«0 /FB , f , k , 2 SIG 

p^nt E, ‘ KEAN p - 

PRINT 
GO TO 110 
END 


S.E 

(P+SE) = ## , ## 



C , 4 TTEST . 

This program is' used to calculate t 


use-.i-as a test -for ' the sicnific. 
means. 


mce of the 


TTKST 


9:96 


FRIDAY AUG 2 4, \ Q7 


nw n-fPUT '■•eam.st ov.no 

I2R IliPUT PS.SS.Ma 
125 INPUT AS 

[30 S3 = ( (Ml -1 )/Nj )*si-T2 
AO M = C(NS-| )/Maj» S2) a 

150 F=Nl+N2-2 

Pf ^“SOFtt <M *S3+US*S4 )/F) 

!asa-i a >» 

190 DEGREES OF F r k v n n m 

200 print* RttDOM 

210 GO TO 110 
220 END 


:he t: 
diffe 
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C . 5 PUT? OL 


This program is used to calculate the 20 t v- 
for each N riin and to store then ih ah’ appropriate * t: 

PUTPOL 9:29 FRIDAY AUG 24, 1973 


i 00 

1 1 0 
120 
J 30 
140 
1 50 
1 60 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
4 m 
410 
420 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
540 


P AMD FILE IT 


REH PUTPOL TO ANALYZE 
UIU 0(300 ),p( 25) 

PKILT OF RUNS TO ANALYZE AND FILE IS 
I NPUT NO 

FOR Nl =| TO N0 
READ I , J,K 
LET C, D=0 
READ 130, T 
FOR M2 =1 TO 80 
READ B 1 ,B2 
C=C+B1/B0 

D=D+B2/B0 1 . 

NEXT M2 
FOR N3 = l TO T 
READ V,N,H,M 
A=(H-C)/(M-D) 

B=(V-C)/<N-D) 

P CN3 ) - ( A-B ) / ( A+B ) 

NEXT N3 

IF J ~ • 2 THEN 330 
PRINT USING 310, J, I 


KEV DATA NOT ADDED TO FILE 


: PRESS=* * • 

GO TO 710 
A$=‘P100* 

IF 1=100 THEN 510 
AS= ' P 150* 

IF 1=150 THEN 510 
AS = ‘ P200 ' 

IF 1=200 THEN 510 
AS = ' P2 50 * 

IF 1=250 THEN 510 
AS = * P300 1 
if 1=300 then 510 

AS= , P350* 

if i =350 then 510 

AS = ' P400 * 

IF I =400 THEN 510 
A S = 1 P 4 5 0 * 

IF 1=450 THEN 510 

PRINT ‘UNACCEPTABLE ENERGY VALUE? 

GO TO 710 

OPEN 1 ,A Sj INPUT 

GET E, G 

IF E=0 THEN 570 
FOR L=1 TO E 


lues 

1 0* • 



55fl GET 0 ( L ) 

560 NEXT L 
570 CLOSE 1 
560 6=1 

590 OPEN 1 * AS / OUTPUT 

600 Q = E+T 

610 PUT 1 

620 FOP. R = 1 TO E 

630 PUT 1 :0(R) 

640 NEXT R 

650 FOR N4=l TO T 

660 PUT 1 :P(N4 ) 

670 NEXT N4 

680 PRINT USING 690 /i 

690 :###'* KEV DATA ADDED TO FILE 

700 CLOSE I 

710 NEXT N I 

720 END 



C. 6 PUT? CLP, 


-•.his program is use ;1 to calculate the 20 


va 1 


for each 3 run. and to store' thee, in a approori ate ' tile 


PUTPOLR 


9;3? FH I DAY AUG 24, 1973 


1 00 
110 
1 20 
130 
140 
150 
160 
170 
1 80 
1 90 
200 
21 0 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
4 J 0 
420 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 

54 0 


IT 


RF» PUTPOL TO ANALYZE P AND FILE 
UIH O(300) / p(25j 

PRINT ■# OF RUNS TO ANALY7E AND FILE 


N0 


B0 


INPUT N0 
FOR Ml = J T(j 
head I, j,k 
L ET C , D = 0 

head B0, T 

i OR N2- 1 TO 
HEAD B l j B2 
C = C-f 81 /B0 
D=D+B2/8C 
NEXT N2 
FOR N3=l TO T 
HEAD v,n,h,m 

CH-C >/( M-D) 
B=(V-C)/CN-D) 

P ( N 3 ) = ( a -B } / ( A + B ) 
NEXT N3 

IF J “ • 2 THEN 330 
PHINT USING 310, J, i 
SPR.ESS si?* • # * 

UO TO 710 
A4«*P100R * 

if i * 1 0 0 then 

AS sr • p j 50 R ' 

IF 1=1 50 THEN 
AS= 'P200R • 

IF 1=200 THEN 
AS = *.P2 50R ' 

IF 1=250 THEN 
AS *=,* P300R ' 

IF 1=300 THEN 
As^'pasoR* 
if 1=350 then 

AS = * P400R » : 

IF 1=400 THEN 
AS = * Pa 50R * 

IF 1=450 THEN 


IS 


KEV DATA NOT ADDED TO FILE 


510 


51 0 
510 
510 
510 


510 


510 


go'tS ;^ accep ^ble energy value- 

OPEN 1, AS, INPUT 
GET E , G 

IF E=0 THEN 570 
FOR L=1 TO E 
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550 GET 0(L) 

560 NEXT L 
570 CLOSE i 
500 G- 1 

590 OPEN 1 , AS , OUTPUT 
600 Q -E+T 
610 PUT I t Q , G 
620 FOR R=i TO E 
630 PUT 1 : 0(R ) 

640 NEXT R 

650 FOR N4 = | TO T 

660 PUT 1 : P <N4 ) 

670 NEXT N4 

680 PRINT USING 690,1 

690 tee#* KEV DATA ADDED TO FILE 

700 CLOSE 1 

710 NEXT N 1 

720 END 



C . 7 .PUTPOLV 


This program is used' to calculate the 20 - vali: 
or each run and to store them according to pressu 

PUTP0LP 9:35 FRIDAY AUG 24, 1 Q73 


IR0 REh PUTPOL to analyse 
i 1 fl DIH O<3<?0),P(2S> 

120 PRINT V* OF RUNS TO 
130 INPUT Up 
140 FOR Nlrsj 7Q 
150 RLAU I , *J , K 


P AN D FILE IT 


ANALYZE AND FILE 


nu 


160 LET C,D = 0 
170 READ B0,T 
180 FOP. N2 = l TO BO 
1 90 READ B l , B2 
200 C=C+B1/B0 
210 D=?D-»-B2/B0 
220 NEXT W2 
230 FOR N3 = i TO T 
240 READ 

250 A = (H-C ) / C M-D ) 

260 B = < V-C >/.(N-D) 

270 PCN3 > = CA-B)/<A+B> 

280 NEXT N3 
290 IF J=.0| THEN 380 
300 IF J= • 0 5 THEN 470 
310 IF J = ♦ 2 THEN 540 
320 IF J = « 6 THEN 570 
330 IF J-1.0 THEN 640 
340 IF’ <J = ] .5 THEN 710 
350 PRINT USING 360, I, J 
360 :PRESS OUT OF BOUNDS 
370 GO TO 980 
380 AS® * T 150 P0 i • 

390 IF 1=150 THEN 780 
400 AS « * T300P0 1 • 

410 IF 1*300 THEN 780 
420 AS = * T 4 50 P0 1 * 

430 IF I =4 50 THEM 7 80 
440 PRINT USING 450, I, J 
450 .-ENERGY OUT OF BOUNDS ###, 
460 GO TO 980 


470 A $ = * T 1 5 0 P 0 5 * 

480 IF 1=150 THEN 780 
490 A S = 1 T 3 0 0 P 0 5 ' 

500 IF 1=300 THEN 780 
510 AS = ' T450P0 5 * 

520 IF I =450 THEN 760 

530 GO TO 440 

540 PRINT USING 550, I, J 





:es 

ener 



130 


55(1 sUSE PUTPOL 000, ##. 
560 GO TO 980. 

570 A $ = * T l 5 0 P 6 * 

560 IF I s 1 50 THEN 780 
590 AS»' T3~POP6 ' 

600 IF 1=300 THEN 780 
610 A$=*T450P6' 

620 IF 1*450 THEN 780 
630 GO TO 440 
640 AS = ' T 1 5-0 P I 0 • 

650 IF 1 = 150 THEN 780 
660 A3> = *T300P1 0' 

670 IF 1=300 THEN 780 
680 A£ = ' T450P 10*. 

690 IF i =450 THEN 7R0 ' 

700 GO TO 440 
7*0 Ai = * T i 50P ] 5 1 
720 IF 1=150 THEN 780 
730 AS= ' T300P | 5 ' 

740 IF 1=300 THEN 780 
750 A£ = *T450P I 5* 

760 IF I =450 THEN 780 

770 GO TO 440 

780 OPEN J, AS, INPUT * 

790 GET E,G 
800 IF E = 0 THEN 840 
810 FOR L=1 TO E 
820 GET 0 ( L ) 

830 NEXT L 
840 CLOSE I 
850 G = 1 

860 OPEN I /AS* OUTPUT 
870 Q =E+T 
880 PUT l:Q,G 
890 FOR R=J TO E 
900 PUT | i 0 (R ) 

910 NEXT R 

920 FOR N4=l TO T 

930 PUT 1 : p < N4 ) 

940 NEXT N4 

950 PRINT USING 960,1,0, AS 
960 t0000 KEV, 00.04 U HE 
970 CLOSE I 
980 NEXT N I 
990 END 


4 0 


data added to 



00000000 




181 


C « 3 PUT 

This program is used to initialize storage fi 
ithin .the computer for- later acceptance of calculate 
. values . 


PUT 


9:39 FRIDAY 


100 PRINT 
J 1 0 INPUT 
120 OPEN \ 
130 M=0 
140 N=I 


'FILE WANE* 
AS> 

OUTPUT 


1^0 PUT- I ,m #N 
160 CLOSE I 
170 GO TO 100 
180 END 


AUG 24, 1973 


les 



102 


C, 9 ? 

xhio program is used to calculate and store :i as 
directed by the user and was used for low current runs. 


PUTPTES T 


9:40 FRIDAY AUG 24, 1973 


DIM 0(300 )/P (25) 

PRINT «# OF RUNS TO ANALY? E AND FILE 
INPUT N0 

FOR N 1 = 1 TO 'm 
READ I,J,K \ 
l»ET C , D~0 
READ B0,T 
FOR W 2 - 1 TO B0 
READ B1,B2 
C=C + Hl /m 
D=D+B2/B0 
NEXT N2 
FOR N3=| TO T 
READ V,N,H.»M 
A=(H-C)/CM-D) 

Ba< v-c )/<N-D) 

P(N3 ) s <A-B) / (A+B) 

N.EXT 


IS • 


c.p; 

.2,90 

h%% 

H V) M 

HI 

K J 0 

33*0 




1# 

mk 

It 

*3 $3 

m 

W 

,420 

V .'Jr* 

« 

i* 

m 

m 

490 

,-600 


< i|^S^‘^ RUN - N0 * ## ,s 
Jfefe V( tea m r 

e 

Wr 

{■> tf! U 

CLuS E„ 1 

. Ztf'J jW 1! ^ AS * OUTPUT 

p^Nj l f AS, OUTPUT 
Q=E+T S 43 

?. '’a" y. 117 

PUT T :/),(, i; ' 

* 


'6k te'PTo 


N F X V "r,T*y 

» 
mw ' 

NitxT Nl 


T 

960 


■AS 


^W/rijfAS 

»*•»* U HE 


•OAT A ADDED TO 

data added to 


$■$ $ A A «? c? :;i 


520 END 



APPENDIX D 


QUALITATIVE CONSIDERATIONS OF THE 
EFFECT OF HELIUM PRESSURE ON tt 

In previous experimental studies of the polariza- 
tion of the light from the 3 1 P-2 1 S transition in Helium, 
the researchers have either claimed that their result was 
independent of the pressure of Helium in the target cham- 
ber below .2 mtorr ( Van Eck (1964), and Van den Bos 
(1969)), or they show a pressure dependence of it down to 
their lowest pressure of .2 mtorr ( Scharmann (1969)). We 
have clearly shown that tt is still pressure dependent 
below .2 mtorr. The problem which remains however is to 
explain that pressure dependence, and especially to ex- 
plain the change in the sign of tt as the He pressure is 
lowered. Scharmann (1969) also observed a sign change, 
but was not able to explain it. 

In order to attempt to discover the reason for the 
sign change in tt , we followed a suggestion made by Pro- 
fessor R. H. Lambert. We reanalyzed our data so as to 
show relative intensities of X 5016 A light as a function 
of He gas pressure. Both the parallel and perpendicular 
intensities were calculated by equation 4.1.2. We then 
normalized these results to He pressure. The results of 
these calculations are shown in Tables Dl, D2 and D3. 



We have plotted the beam current and pressure- normalized 
results in Figures Dl, D2 and D3. These intensities are 
proportional to the number of photons per proton per 
target atom received at the detector. 

We had intended to look into changes in the rela- 
tive intensities of parallel and perpendicularly polarized 
light (as a function of pressure) in order to attempt an 
explanation of the change in sign of tt . However the 
relative changes in relative intensities are very small 
compared to the average relative intensity changes (as a 
function of pressure)., therefore we were not able to find 
a probable cause for the change in sign of tt. 

The information gained however, has caused us to 
realize that significant trapping in the metastable states 
of He occurs within the target chamber. We have pre- 
viously discussed (section 4.2) trapping at the level 
and how it would tend to depolarize, the A 5016 A light. 

We have also previously corrected for this effect. How- 
ever capture at levels other than the ground state must- 
now be considered. Inspection of figure 11 . shows three 
states which do not connect (via electric dipole transi- 
tions) with the 1 J S ground state. These three states 
(which are long lived and called metastable) are the 2 1 s, 

2 3 S an <* the 2 3 p. The triplet states need not concern us 
since transitions between triplet states will not result 
in light of wavelength 5016 A, and since the only inter- 
system transition of note is 2 3 P-l^s which simply leaves 



TABLE Dl 

INTENSITY Vs He PRESSURE AT 150 keV BEAM ENERGY 

NORMALIZED TO BOTH 

AV. NORMALIZED TO BEAM CURRENT AND 


He 

PRESSURE 

mtorr 

BEAM 

CURRENT 

UA 

# OF 
RUNS 

# OF 

MEASURE- 

MENTS 

S.E. 

OF 

1// 

BEAM 

MEAN 

X,/ 

CURRENT 

MEAN 

1^ 

S.E. 

OF 

He 

MEAN 

I// 

PRESSURE 

MEAN 

.01 

5.5 

8 

160 

.0015 

.082 

.083 

.0014 

8.2 ±.15 

8.3 ±.14 

.05 

6.6 

5 

100 

.011 

.276 

.261 

.01 

5 . 52± . 22 

5 . 2 2 ± . 2 

.2 

7.1 

7 

140 

.012 

.933 

.898 

.011 

4 . 66 ± . 06 

4 . 48 ±. 06 

.6 

7. 

2 

40 

.051 

4.133 

4.082 

.046 

6 . 88 ± . 08 

6.8 ±.08 

1.0 

5.8 

2 

40 

. .037 

8.770 

8.706 

.026 

8 . 77 ± . 04 

8.71±.03 

1.5 

6. 

2 

40 

.153 

16.759 

16.99 

.161 

11.0 3 ±.10 

11. 3 1 ± . 1 1 
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TABLE D2 


INTENSITY Vs He PRESSURE AT 300 keV BEAM ENERGY 


He 

PRESSURE 

mtorr 

AV. 

BEAM 

CURRENT 

UA 

# OF 
RUNS 

# OF 

MEASURE- 

MENTS 

S.E. 

OF 

• I// 

.01 

11.8 

7 

140 

.001 

.05 

11 . 4 

7 

140 

.002 

.2 

10.5 

10 

200 

.004 

.6 

12.2 

3 

60 

.106 

1.0 

11.2 

3 

60 

.323 

1.5 

11.3 

2 

40 

.662 


NORMALIZED TO BOTH 

NORMALIZED TO BEAM CURRENT AND 


BEAM 

MEAN 

V 

CURRENT 

MEAN 

S.E. 

OF 

I_L 

He 

MEAN 

I, 

PRESSURE 

MEAN 

.0629 


.0615 

.0009 

6.291.1 

6.151.09 

.159 


.162 

.002 

3.191.04 

3.241.04 

.594 


.624 

.004 

2.971.02 

3.121.02 

1.907 

2 

.06 

.118 

3.181.18 

3.441.2 

4.428 

4 

.714 

.353 

4.431.32 

4.711.35 

11.917 

12 

.528 

.717 

7.951.44 

8.381.48 
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TABLE D3 


INTENSITY Vs He PRESSURE AT 450 keV BEAM ENERGY 


NORMALIZED TO BOTH 

NORMALIZED TO BEAM CURRENT AND 


He 

PRESSURE 

mtorr 

* * ¥ • 
BEAM 
CURRENT 
UA 

# OF 
RUNS 

# OF 

MEASURE- 

MENTS 

S.E. 

OF 

*// 

BEAM CURRENT 
MEAN MEAN 

Z // 

S.E. 

OF 

. ; , He 

MEAN 

: '' 1 // 

PRESSURE 

MEAN 

IjL 

.01 

11.8 

4 

80 

.0018 

.0877 

.0914 

.0018 

6.771.18 

9.141.18 

.05 

10.4 

5 

100 

.0029 

.1741 

.1813 

.0023 

3.481.58 

3.621.46 

.2 

10.3 

7 

140 

.004 

.4821 

. 5387 

.005 

2.411.02 

2.691.03 

.6 

10.7 

3 

60 

.039 

1.953 

2.181 

.058 

3.261,07 

3.641.1 

1.0 

10.8 

3 

60 

.171 

3.754 

4.094 

.211 

3.751.17 

4.091.21 

1.5 

10.5 

2 

40 

.106 

5.107 

5.47 

.129 

3.411.07 

3.651.09 
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INTENSITY Vs HELIUM PRESSURE 
@ 150 keV BEAM ENERGY 



0 .2 .6 1.0 1.5 

He Pressure mtorr 


o= Z„ ' lx 


FIGURE Dl 
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INTENSITY Vs HELIUM PRESSURE 
<5> 300 keV BEAM ENERGY 
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FIGURE D2 
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INTENSITY Vs HELIUM PRESSURE 
@ 450 keV BEAM ENERGY 



He Pressure mtorr 

« * l l 


FIGURE D3 
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a ground state atom. It is then the 2 l s state with which 
we must concern ourselves. In order to estimate the 
effect on A 5016 A photons by 2 ! S state He atoms, we must 
first discuss the presence of 2 ! S atoms in the target 
chamber. Referring to calculations done in sections 4.2 
and 4.3, we have at our highest working pressure (1.5 
mtorr) a mean free path of 47 cm, a mean collision time 
of ^.3 msec, and an estimated minimum wall collision time 
of ^.02 msec. Similar calculations for our lowest work- 
ing pressure (.01 mtorr) yield a mean free path of 
000 cm, a mean collision time of v50 msec, and an 
estimated minimum wall collision time of a,, 02 msec. 

Note that the wall collision times are minimum times and 
in fact the "diffusion" times would be much longer. 
Furthermore, we have from Hasted (197.2) the fact that the 
natural lifetime of the metastable state He 2 1 S is 
> 1 msec. We can then conclude that at our high pressure 
runs the He 2 ! S level will be collisionally quenched and 
significant absorbtion of A 5016 A light will not occur. 
However as we go to lower pressures, the number of sur- 
viving He 2 1 S atoms increases and we can expect some 
absorbtion of the A 5016 A light. Finally at very low 
pressures one would expect a high percentage of He 2 l S 
in the target chamber, but these metastables would be 
spread so thinly that there would be lit e photon- 
metastable interaction, and therefore ar. increase in the 
amount of A 5016 A light detected. 
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The results which have been normalized to both 
Beam Current and Gas Pressure are the relative intensi- 
ties of A 5016 A light per proton per target atom. The 
minimum observed in these curves is just what one would 
expect (qualitatively) from our discussion of the exis- 
tence of the metastable He 2 l S atoms. We thus conclude 
that trapping by the metastable He 2 1 S atoms is respon- 
sible for the changes in relative intensity with pres- 
sure. 

In conclusion then, we are still unable to ex- 
plain the cause of the change in the sign of it as the 
He pressure changes. We have shown qualitatively that 
trapping is occurring at the metastable 2 X S level. This 
trapping provides yet another reason for the next experi- 
menter to work at lower pressures. Only then will the 
photon from a beam excited 3 1 P-2 1 S transition be directly 


detected . 



